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Performance  of  Nonlinear  Mechanical,  Resonant-Shunted  Piezoelectric,  and 
Electronic  Vibration  Absorbers  for  Multi-Degree-of-Freedom  Structures 


Gregory  S.  Agnes 


(ABSTRACT) 


Linear  vibration  absorbers  are  a  valuable  tool  used  to  suppress  vibrations  due  to  har¬ 
monic  excitation  in  structural  systems.  Limited  evaluation  of  the  performance  of  nonlinear 
vibration  absorbers  for  nonlinear  structures  exists  in  the  current  literature.  The  state  of 
the  art  is  extended  in  this  work  to  vibration  absorbers  in  their  three  major  physical  imple¬ 
mentations:  the  mechanical  vibration  absorber,  the  inductive-resistive  shunted  piezoelectric 
vibration  absorber,  and  the  electronic  vibration  absorber  (also  denoted  a  positive  position 
feedback  controller).  A  single,  consistent,  physically  similar  model  capable  of  examining 
the  response  of  all  three  devices  is  developed. 

The  performance  of  vibration  absorbers  attached  to  single-degree-of-freedom  structures 
is  next  examined  for  performance,  robustness,  and  stability.  Perturbation  techniques  and 
numerical  analysis  combine  to  yield  insight  into  the  tuning  of  nonlinear  vibration  absorbers 
for  both  linear  and  nonlinear  structures.  The  results  both  clarify  and  validate  the  existing 
literature  on  mechanical  vibration  absorbers.  Several  new  results,  including  an  analytical 
expression  for  the  suppression  region’s  location  and  bandwidth  and  requirements  for  its 
robust  performance,  are  derived. 

Nonlinear  multiple-degree-of-freedom  structures  are  next  evaluated.  The  theory  of  Non¬ 
linear  Normal  Modes  is  extended  to  include  consideration  of  modal  damping,  excitation, 
and  small  linear  coupling,  allowing  estimation  of  vibration  absorber  performance.  The  dy¬ 
namics  of  the  N-|-l-degree-of- freedom  system  reduce  to  those  of  a  two-degree-of-freedom 
system  on  a  four-dimensional  nonlinear  modal  manifold,  thereby  simplifying  the  analysis. 
Quantitative  agreement  is  shown  to  require  a  higher  order  model  which  is  recommended  for 
future  investigation. 

Finally,  experimental  investigation  on  both  single  and  multi-degree-of-freedom  systems 
is  performed  since  few  experiments  on  this  topic  are  reported  in  the  literature.  The  ex¬ 
perimental  results  qualitatively  verify  the  analytical  models  derived  in  this  work.  The 
dissertation  concludes  with  a  discussion  of  future  work  which  remains  to  allow  nonlinear 
vibration  absorbers,  in  all  three  physical  implementations,  to  enter  the  engineer’s  toolbox. 
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Chapter  1 


Introduction 


Suppression  of  vibration  is  increasingly  demanded  by  modern  structures.  The  drive  to 
produce  higher  performance,  lower  cost  structures  and  machines  has  resulted  in  reduced 
structural  weight  and  hence  increased  compliance.  Nowhere  has  this  trend  been  more 
evident  than  the  aerospace  industry  where  the  demand  for  smaller,  lighter,  cheaper  aircraft 
and  spacecraft  systems  drives  the  acquisition  effort.  The  demand  for  performance  exists 
in  the  face  of  ever  harsher  environments  with  increased  disturbance  levels.  Meanwhile,  the 
need  for  cleaner,  more  maintainable  systems  reduces  the  number  of  connectors,  joints,  wires, 
and  other  “dirty”  components  which  contribute  to  structural  damping.  Modern  systems 
thus  face  higher  disturbances  with  lower  damping  and  increased  compliance,  resulting  in 
larger  and  often  nonlinear  vibrations.  To  achieve  performance  goals,  the  designer  must 
address  the  suppression  of  nonlinear  vibrations. 

1.1  Review  of  the  Vibration  Suppression  Literature 

Techniques  used  to  reduce  structural  vibrations  can  be  classified  with  two  criteria.  First, 
one  can  either  modify  existing  modes  or  add  additional  modes  to  the  system.  Secondly,  the 
technique  can  be  either  active  or  passive.  A  technique  is  passive  if  it  requires  no  external 
energy;  otherwise,  it  is  active.  Differentiating  between  adding  modes  or  altering  existing 


1 


modes  results  in  an  important  distinction  in  the  response  of  the  modified  system.  This  is 
not  well  recognized  in  the  literature  or  in  practice  and  can  lead  to  valuable  insight  into  the 
design  and  implementation  of  vibration  suppression  systems.  The  purpose  of  this  section  is 
to  review  the  relevant  literature  within  this  framework,  motivating  the  subsequent  research. 

1.1.1  Modal  Modification 

Passive  modal  modification  is  perhaps  the  simplest  and  most  commonly  used  technique. 
To  reduce  the  vibrations,  the  mode  shapes  and  modal  parameters  are  altered  until  accept¬ 
able  performance  is  achieved.  For  instance,  if  the  rotation  frequency  of  machinery  is  known 
to  lie  within  a  certain  bandwidth,  the  support  structure  is  designed  such  that  no  resonant 
frequencies  lie  within  that  bandwidth.  Another  example  of  passive  modal  modification  is 
the  use  of  passive  damping  treatments  such  as  constrained  layer  damping  treatments  [1] 
or  the  addition  of  dashpot  dampers  to  a  structure.  These  modify  the  modal  damping  of 
various  structural  modes  and  reduce  the  vibration  of  the  system.  The  performance  of  pas¬ 
sive  modal  modification  can  be  assessed  for  linear  systems  using  techniques  found  in  most 
vibration  texts  [2,3].  First,  the  governing  equations  are  transformed  using  modal  analysis 
into  a  set  of  uncoupled  single  degree  of  freedom  modal  equations,  assuming  proportional 
damping.  The  response  of  each  mode  is  then  determined  and  compared  to  design  criteria. 
Structural  modifications  can  then  be  designed  and  the  structure  reanalyzed. 

Nonlinearity  complicates  this  approach  since  structural  modes  can  respond  to  excitation 
far  away  from  their  resonant  frequencies.  For  example: 

“von  Karman  observed  that  certain  parts  of  an  airplane  can  be  violently  excited 
by  an  engine  running  at  an  angular  speed  much  larger  than  their  natural  fre¬ 
quencies,  and  Lefschetz  described  a  commercial  aircraft  in  which  the  propellers 
induced  a  subharmonic  vibration  in  the  wings  which  in  turn  induced  a  subhar¬ 
monic  vibration  in  the  rudder.  The  oscillations  were  violent  enough  to  cause 
tragic  consequences.”  [4] 
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Nonlinear  systems  can  also  exhibit  jump  phenomena  and  either  quasiperiodic  or  chaotic 
responses  to  harmonic  disturbance.  The  effect  of  changes  in  various  parameters  can  be 
examined  to  a  limited  extent  using  perturbation  techniques  for  either  weak  nonlinearities 
or  small  motions.  Nayfeh  and  Mook  [4]  discuss  such  techniques,  but  their  applications  are 
mostly  limited  to  single-degree-of-freedom  systems.  Multiple-degree-of-£reedom  structures 
axe  treated,  but  not  in  a  modal  manner,  thus  limiting  the  analysis  to  systems  with  a  limited 
number  of  degrees  of  freedom. 

The  concept  of  modal  behavior  for  nonlinear  systems  was  introduced  by  Rosenberg  [5- 
7].  According  to  his  definition,  a  nonlinear  normal  mode  of  a  discrete  conservative  system 
must  satisfy  three  conditions  [8,9]: 

i.  The  masses  execute  periodic  motions  (not  necessarily  harmonic). 

ii.  The  masses  pass  through  static  equilibrium  simultaneously. 

iii.  Displacement  of  one  mass  uniquely  defines  the  displacements  of  all  masses. 

These  criteria  are  analogous  to  the  synchronous  motion  criteria  of  linear  normal  modes. 
However,  the  nonlinear  normal  modes  are  also  classified  as  similar  (the  relation  of  criterion 
three  is  linear)  or  non-similar  (the  relation  is  nonlinear)  in  contrast  to  linear  normal  modes 
which  always  have  linear  relationships  represented  by  the  mode  shape.  The  number  of 
nonlinear  normal  modes  of  a  system  under  this  definition  is  at  least  equal  to  the  number  of 
degrees  of  freedom,  but  can  be  greater. 

Shaw  and  Pierre  [10, 11]  present  an  alternative  definition  of  nonlinear  normal  modes, 
defined  for  both  conservative  and  dissipative  systems.  Their  definition  is  made  in  the  phase 
space  and  the  nonlinear  normal  modes  are  as  follows  [8]: 

i.  Motions  (not  necessarily  periodic)  exist  on  invariant  manifolds  of  the  system  which 
are  tangent  to  the  linear  normal  modes  at  the  equilibrium  points. 

ii.  Motions  pass  through  stable  equilibrium. 
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iii.  Displacement  of  one  mass  uniquely  defines  the  displacements  of  all  masses. 

Thus  the  first  criterion  ensures  there  are  exactly  as  many  nonlinear  normal  modes  as  there 
are  linear  normal  modes  for  infinitesimal  nonlinearity  or  system  energy.  Again,  both  similar 
and  nonsimilar  modes  are  possible.  This  concept  was  later  broadened  by  Slater  and  Inman 
[12]  who  removed  the  second  constraint.  This  allows  nonlinear  normal  modes  to  be  the 
counterpart  of  linear  complex  modes  for  systems  which  are  nonmodally  damped.  Such 
modes  are  designated  nonequal  phase.  Using  the  invariant  manifold  concept,  nonlinear 
normal  modes  can  be  used  to  define  a  nonlinear  coordinate  transformation  from  the  physical 
to  the  modal  system  [11, 12].  In  the  modal  space,  however,  nonlinear  couplings  between  the 
modes  can  exist,  but  in  the  absence  of  internal  resonance,  are  higher  order  eflFects  near  modal 
resonance.  Thus,  nonlinear  normal  modes  can  be  used  to  examine  the  modal  response  of  a 
system,  and  the  results  can  be  used  to  design  passive  modal  modifications  to  the  structure 
in  a  manner  similar  to  linear  normal  modes. 

The  solution  for  nonlinear  normal  modes  (as  defined  by  Shaw  and  Pierre)  can  be  per¬ 
formed  using  a  number  of  techniques  [8].  The  original  Real- Variable  Invariant-Manifold 
method  [10],  while  conceptually  significant,  is  difiicult  to  implement  for  large  systems,  es¬ 
pecially  in  the  presence  of  internal  resonance.  It  can  be  simplified  by  solving  the  linear 
eigenvalue  problem  and  transforming  the  equations  to  modal  space.  Two  methods  pro¬ 
posed  in  [13]  are  easier  to  implement:  the  Complex- Variable  Invariant-Manifold  approach 
and  the  Method  of  Multiple  Scales.  Solution  of  the  manifold  dynamics  for  any  of  these 
methods  is  limited  to  small  nonlinearities  except  for  special  cases. 

Active  modal  modification  usually  involves  feedback  control  using  measurements  at  one 
or  more  points  on  the  structure.  Through  feedback  control,  the  modal  parameters  of  the 
closed  loop  system  can  be  altered.  For  linear  systems,  the  use  of  collocated  velocity  feed¬ 
back  to  provide  increased  structural  damping  is  an  example.  For  multi-degree-of-freedom 
systems,  any  static  gain  feedback  methodology  such  as  eigenstructure  assignment  or  LQG 
control  also  simply  modifies  the  eigenstructure  (modal  structure)  of  the  system.  Many 
powerful  techniques  have  been  developed  for  linear  systems  control  design.  One  of  the  sim- 
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plest  theoretical  methods  is  modal  control  which  determines  changes  in  the  modal  dynamics 
required  and  then  uses  the  inverse  modal  transform  to  determine  the  physical  controller. 
Perfect  modal  observability  and  controllability  are  required  to  implement  this  approach. 

For  nonlinear  systems,  Slater  [12, 14]  extended  the  approach  of  Shaw  and  Pierre  [11]  to 
include  forcing  and  sensing.  Thus  nonlinear  controllers  could  be  designed  in  the  nonlinear 
modal  coordinate  space  and  then  transformed  into  physical  coordinates.  In  a  similar  man¬ 
ner,  individual  modes  could  be  linearized,  creating  a  modal  form  of  feedback  linearization 
control.  In  contrast  to  linear  theory,  nonlinear  theory  lacks  the  tools  for  optimal  design, 
requiring  a  trial  and  error  approach.  Whether  active  or  passive  modal  modification  is 
employed,  the  result  is  a  physical  structure  with  an  improved  modal  structure  capable  of 
satisfying  performance  requirements. 

1.1.2  Modal  Addition 

When  modal  modification  is  impractical,  additional  degrees  of  freedom  can  be  placed 
into  the  structure  in  such  a  manner  that  vibration  of  critical  components  is  reduced.  Usu¬ 
ally  these  new  degrees  of  freedom  will  undergo  large  motions,  absorbing  the  energy  of  the 
disturbance.  For  example,  vibration  isolation  is  often  achieved  by  placing  a  soft  spring  in 
the  energy  transmission  path.  If  a  hard  mount  is  replaced  by  a  soft  spring,  a  new  degree  of 
freedom  is  added  to  the  structure.  (If  an  existing  flexible  mount  is  modified,  this  would  be 
modal  modification  instead.)  This  spring  will  undergo  large  deflection,  but  isolates  the  crit¬ 
ical  system  components  from  the  disturbance.  Both  passive  and  active  vibration  isolation 
are  possible  (see  Sciulli  [15]  for  a  recent  review  of  vibration  isolation  for  both  single  and 
multiple  degree  of  freedom  systems).  Nonlinear  effects  can  be  used  to  improve  vibration 
isolation  particularly  the  concept  of  mode  localization  [9]. 

When  it  is  not  possible  to  isolate  a  system  from  disturbances,  vibration  absorbers  use 
modal  addition  to  passively  attenuate  the  response  of  a  structure  to  external  excitations. 
The  mechanical  vibration  absorber  Weis  first  introduced  at  the  turn  of  the  century.  The 
theory  of  the  classical  device  may  be  found  in  many  standard  vibrations  textbooks.  See  for 
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instance  den  Hartog  [16]  or  Inman  [2].  In  brief,  a  single-degree-of-freedom  system  is  inertially 
coupled  to  a  structure  with  its  resonance  tuned  near  that  of  the  primary  structure,  creating 
a  one-to-one  frequency  relationship.  Properly  tuned,  the  device  can  reduce  the  response  of 
a  structure  to  a  narrowband  harmonic  excitation  or,  by  adding  damping  to  the  absorber, 
to  a  broadband  excitation. 

A  recent  review  of  the  theory  and  application  of  the  linear  vibration  absorber  [17]  notes 
the  mass  and  stroke  length  design  limitations  imposed  by  practical  applications.  Geometric 
considerations  can  limit  the  stroke  length  (the  amount  of  travel)  of  the  added  mass.  Since 
the  mciss  added  to  a  system  must  often  be  minimized,  the  tradeoff  between  mass  and  stroke  is 
a  crucial  design  constraint.  Given  a  fixed  mass  and  stroke  length,  an  absorber  is  efficient  for 
only  a  narrow  range  of  driving  force  and  fi'equency  variations.  Consequently,  active  and/or 
adaptive  feedback  control  is  topical  [17, 18].  Using  active  control  to  minimize  the  response 
of  the  system,  smaller  absorber  stroke  lengths  can  be  achieved,  albeit  with  increased  system 
complexity. 

One  solution  to  these  limitations  is  the  piezoelectric  vibration  absorber  formed  by  cou¬ 
pling  a  resonant  electrical  shunt  between  the  electrodes  of  piezoelectric  materials  attached 
or  embedded  in  the  structure.  Tuning  the  shunt  near  a  structural  mode  causes  mechanical 
energy  to  be  transformed  into  electrical  energy  and  dissipated  by  the  resistive  element  in 
the  shunt  in  a  manner  analogous  to  a  damped  mechanical  vibration  absorber.  Forward 
[19]  experimentally  demonstrated  the  use  of  both  resistive  and  resonant  electrical  shunts. 
Hagood  and  Von  Flotow  [20]  later  presented  an  analytical  model.  Other  researchers  [21-23] 
have  extended  this  device  providing  multimode  and  self-  tuning  shunts.  Incorporation  of  a 
simultaneous  active  and  pa.ssive  signal  is  also  possible  [24]. 

The  piezoelectric  vibration  absorber  solves  the  stroke  length  problem  by  providing  an 
unobtrusive  energy  dissipation  mechanism,  but  creates  new  design  challenges.  The  large 
inductances  required  to  tune  the  electrical  resonance  axe  unreasonable  for  passive  compo¬ 
nents.  Consequently,  active  inductors  are  synthesized  using  op-  amps  making  piezoelectric 
vibration  absorbers  semiactive  (or  powered)  devices.  A  second  design  challenge  presented 
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by  the  piezoelectric  vibration  absorber  is  the  structural  properties  of  available  piezoelectric 
materials.  Since  piezoelectric  materials  transform  strain  energy  into  electrical  energy,  they 
must  be  located  in  areas  of  high  strain  energy  to  function  efficiently.  Since  their  structural 
properties  are  inferior  to  traditional  load  carrying  materials,  this  poses  a  design  challenge: 
maximize  strain  energy  in  the  piezoelectric  material  while  maintaining  structural  integrity. 
A  final  concern  for  the  designers  is  the  nonlinear  behavior  of  piezoceramics  [25].  To  date,  no 
study  could  be  found  on  the  effects  of  small  nonlinearities  on  the  performance  of  piezoelectric 
vibration  absorbers. 

Results  have,  however,  been  reported  on  nonlinear  (mechanical)  vibration  absorbers. 
Early  work  [26-28]  considered  a  linear  system  with  a  nonlinear  absorber  and  neglected 
damping  entirely.  Most  commonly,  a  nonlinear  spring  is  used  in  place  of  the  linear  spring 
when  designing  the  absorber.  In  1952,  Roberson  [26]  reported  that  a  broadening  of  the 
suppression  bandwidth  was  possible,  and  that  the  zero  response  frequency  (design  opera¬ 
tional  frequency  in  most  cases)  changed  with  force  level  and  nonlinearity.  References  [27, 28] 
came  to  similar  conclusions  with  slightly  different  systems.  Miller  [29]  presented  a  nonlinear 
design  using  pneumatic  springs,  while  Hunt  and  Nissen  [30]  used  Belleville  washers  and  re¬ 
ported  that  a  broadening  of  the  suppression  bandwidth  is  possible.  Several  authors  [31-33] 
have  applied  optimization  methods  to  absorber  design,  but  considered  system  damping  of 
ten  percent  or  more.  This  is  unrealistic  for  many  applications.  Each  pointed  to  bandwidth 
improvements  and  several  claimed  higher  harmonic  behavior  improved  the  design. 

However,  Shaw  [34]  warned  of  the  presence  of  quasiperiodic  responses  near  the  tuning 
frequency.  Such  motions  were  present  in  numerical  simulations  of  lightly  damped  linear 
systems  with  nonlinear  absorbers.  These  responses  were  also  confirmed  analytically  in 
[35, 36].  Rice  and  McCraith  [35]  used  the  method  of  harmonic  balance  to  numerically  inves¬ 
tigate  the  response  of  a  linear  system  with  a  nonlinear  absorber  to  narrowband  excitation. 
Natsiavas  [36]  considered  the  case  of  the  nonlinear  system  with  a  nonlinear  absorber  and 
demonstrated  the  adverse  effects  of  system  nonlinearity  on  absorber  performance.  Note 
that  limited  experimental  work  has  been  reported. 
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Alternately,  the  motion  of  the  structure  can  serve  as  a  parametric  excitation  to  a  non¬ 
linear  absorber  [37].  Oueini,  et  al.  [38]  have  experimentally  implemented  a  similar  system 
using  nonlinear  feedback  control  with  piezoelectric  actuators.  This  approach  will  not  be 
considered  in  this  work.  Thus  there  appeared  to  be  both  some  promise  and  some  peril  in 
using  nonlinear  mechanical  vibration  absorbers.  The  effect  of  nonlinearity  on  the  response 
of  the  piezoelectric  vibration  absorber  is  cmrrently  unknown. 

Active  modal  modification  has  also  been  pursued.  For  a  linear  system,  any  dynamic 
compensator  designed  via  Hqo  ,  p-synthesis,  or  other  modern  control  algorithms  uses  a 
combination  of  modal  modification  and  modal  addition,  but  physical  intuition  into  the 
inner  workings  of  the  control  schemes  is  limited  since  they  are  designed  for  generic  first 
order  systems  of  the  form  x  =  Ax  instead  of  second  order  dynamic  systems. 

In  contrast,  the  Positive  Position  Feedback  (PPF)  algorithm  developed  by  Goh  and 
Caughey  [39]  and  implemented  by  Fanson  and  Caughey  [40, 41]  uses  the  second  order  form, 
allowing  better  physical  insight  to  vibration  control  by  active  modal  addition.  In  this  algo¬ 
rithm,  a  position  signal  is  compensated  by  a  second  order  filter  for  feedback  control.  For 
linear  systems,  the  PPF  controller  is  stable  even  in  the  presence  of  unmodelled  actuator 
dynamics,  unlike  the  often  used  direct  velocity  feedback  (which  is  a  form  of  modal  modi¬ 
fication).  Baz,  et  al.,  later  combined  PPF  with  Independent  Modal  Space  Control  [42,43] 
to  design  independent  second  order  feedback  compensators  for  individual  modes.  Caughey 
later  noted  that  positive  position  feedback  was  a  generalization  of  the  mechanical  vibra¬ 
tion  absorber  [44].  For  this  reason,  PPF  control  can  be  denoted  an  electronic  vibration 
absorber.  Other  significant  contributions  to  the  PPF  literature  include  the  work  of  Dosch 
[45],  who  derived  stability  conditions  and  applied  PPF  to  an  eight-ribbed  space  antenna. 
Many  other  numerical  and  experimental  implementations  of  the  PPF  control  scheme  may 
be  found  in  the  literature.  No  reports  on  the  influence  of  nonlinearity  on  the  performance 
of  PPF  controllers  could  be  found. 
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1.2  Objectives 


The  focus  of  this  work  will  be  on  vibration  suppression  for  nonlinear  multi-degree-of- 
freedom  discrete  structures  via  vibration  absorbers.  The  detailed  objectives  are  as  follows: 

•  Demonstrate  the  equivalence,  to  first-order  nonlinear  effects,  of  the  mechanical  vibra¬ 
tion  absorber,  the  piezoelectric  vibration  absorber,  and  positive  position  feedback. 

•  Utilize  perturbation  methods  to  find  the  general  nonlinear  response  of  the  aforemen¬ 
tioned  systems. 

•  Develop  design  guidelines  for  tuning  the  linear  and  nonlinear  parameters  for  vibration 
absorbers  for  nonlinear  single-degree-of-freedom  systems. 

•  Demonstrate  the  use  of  nonlinear  normal  modes  to  allow  modal  design  of  absorbers 
for  N-degree-of-freedom  nonlinear  structures  using  these  guidelines. 

•  Qualitatively  verify  the  models  through  experimental  investigations  for  both  single 
and  multiple  degree-of-freedom  systems. 

Analytical,  numerical,  and  experimental  techniques  will  be  used  to  demonstrate  the  perfor¬ 
mance  of  nonlinear  absorbers  for  a  variety  of  nonlinear  systems. 

1.3  Overview 

To  achieve  the  objectives  presented  in  the  previous  section,  nonlinear  models  for  the 
mechanical,  shunted  piezoelectric  and  electronic  (PPF)  vibration  absorbers  are  derived  in 
Chapter  2.  By  proper  choice  of  states  and  nondimensional  parameters,  these  systems  all 
reduce  to  a  single  set  of  identical  nonlinear  differential  equations  (to  first  order  nonlinear 
effects)  while  maintaining  consistent  physical  variables. 

Chapter  3  is  devoted  to  examining  both  the  linear  and  nonlinear  response  of  the  unified 
model.  Since  there  are  many  control  parameters,  a  limited  set  of  cases  will  be  investi¬ 
gated:  linear  system/linear  absorber  (LSLA),  nonlinear  (cubic)  system/  linear  absorber 
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(NSLA)  and  linear  system/nonlinear  (cubic)  absorber  (LSNA).  Prom  the  work  cited  above, 
quasiperiodic  responses  can  result  from  the  nonlinearities  in  the  LSNA  system.  It  will  be 
shown  analytically  that  they  can  also  occur  in  the  NSLA.  Finally,  a  fourth  system,  non¬ 
linear  (cubic)  system/nonlineax  (cubic)  absorber  (NSNA),  will  be  investigated.  Tuning  is 
accomplished  by  observing  the  effects  of  linear  and  nonlinear  parameters  from  the  previous 
three  systems.  Criteria  for  judging  controller  performance  will  also  be  addressed.  Both 
numerical  and  analytical  (perturbation)  methods  will  be  used  to  calculate  the  responses. 

Once  a  design  methodology  for  a  single-degree-of-freedom  structure  has  been  established, 
N-degree-of-freedom  discrete  systems  will  be  discussed  in  Chapter  4.  The  equations  of 
motion  will  be  placed  into  modal  form  using  nonlinear  normal  modes,  and  an  absorber 
designed  for  one  of  the  modes.  The  procediure  followed  will  be  similar  to  Slater’s,  however 
dynamic  compensation  instead  of  static  compensation  will  be  used.  Since  the  example  of 
Shaw  and  Pierre  [10]  of  a  two-degree-of-freedom  system  has  been  used  often  in  the  literature 
as  a  benchmark  problem,  it  will  be  used  to  illustrate  the  methodology  for  absorber  design. 

Experimental  investigations  will  be  used  to  verify  the  analytical  and  numerical  analysis 
of  the  response  of  the  single-degree-of-freedom  system  with  a  vibration  absorber  using  an 
analog  computer.  These  will  be  conducted  for  various  combinations  of  parameters.  Then 
the  two-degree-of-freedom  system  of  Shaw  and  Pierre  will  be  similarly  investigated.  Finally, 
a  vibration  absorber  will  be  implemented  using  nonlinear  PPF  on  a  beam  to  further  verify 
the  utility  of  the  models  in  an  experimental  application. 

The  proposed  work  will  meet  the  objectives  outlined  above.  A  single  nonlinear  model 
for  a  single-degree-of-freedom  mechanical,  piezoelectric,  or  electronic  vibration  absorber  will 
be  found.  This  model  will  be  used  to  verify  past  research  efforts  for  LSNA  systems  as  well 
as  providing  novel  investigation  into  NSLA  and  NSNA  systems.  The  model  will  preserve 
physical  intuition,  unlike  some  past  investigations.  The  discussion  of  nonlinear  effects  for 
the  piezoelectric  vibration  absorbers  is  new  to  the  literature  and  was  recently  published 
[46].  The  derivation  of  a  tuning  strategy  for  the  design  of  nonlinear  absorbers  will  overcome 
shortcomings  noted  by  previous  research  and  provide  some  analytic  results  for  design.  The 
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state  of  the  art  will  be  further  advanced  by  the  use  of  nonlinear  normal  modes  to  design 
vibration  absorbers  for  N-dof  systems.  Finally,  experimental  verification  of  the  predicted 
quasiperiodic  response  will  be  accomplished.  In  summation,  using  nonlinear  normal  modes, 
a  design  methodology  for  second-order,  nonlinear  vibration  absorbers  will  be  developed  to 
suppress  the  vibration  of  multiple-  degree-of-freedom  structures  and  machines  which  can  be 
implemented  in  mechanical,  shunted  piezoelectric,  or  electronic  (PPF)  form. 
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Chapter  2 


Single-Degree-of-Preedom  System 
Models 


System  models  for  the  mechanical  vibration  absorber,  the  piezoelectric  vibration  ab¬ 
sorber,  and  the  positive  position  feedback  algorithm  may  be  found  in  the  various  references 
cited.  However,  the  common  form  of  the  models  results  in  three  different  sets  of  equations. 
By  careful  selection  of  nondimensional  parameters,  and  by  enforcing  the  requirement  that 
the  degrees  of  freedom  represent  the  motion  of  the  base  system  and  the  stroke  length  of  the 
absorber,  a  common  set  of  equations  can  be  formed. 


2.1  Mechanical  Vibration  Absorber 

A  model  of  a  single-degree-of-freedom  structure  with  a  mechanical  vibration  absorber 
is  shown  in  Figure  2.1.  A  single-degree-of-freedom  system  of  mass,  M,  stiffness,  AT,  and 
viscous  damping,  C,  is  driven  by  an  external  excitation,  F,  and  responds  with  displacement, 
X.  By  adding  an  absorber  (with  mass,  stiffness,  damping  and  displacement  denoted  Mq, 
Ka,  Ca  and  Xa,  respectively)  the  displacement  of  the  original  system  can  be  reduced. 
The  displacement  of  the  absorber  is  also  referred  to  cis  the  stroke  length  in  the  literature. 
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Figure  2.1:  The  mechanical  vibration  absorber. 


Nonlinear  stiffness  of  both  the  system  and  absorber  will  be  considered  of  the  form: 


K{X)  =  KX  +  TX^  (2.1a) 

Ka{Xa)  =  KaXa  +  TaXl  (2.1b) 

Prom  Newton’s  law,  the  equations  of  motion  for  the  mechanical  vibration  absorber  are  then: 

MX  +  CX-  CaXa  ■\-KX-  KaXa  =F  -  TX^  +  TaX^  (2.2a) 

MaX  +  MaXa  +  +  KaXa  =  -  TaX^  (2.2b) 


where  the  overdot  denotes  total  differentiation  with  respect  to  time,  T.  The  term  X  is  then 
solved  for  in  Equation  (2.2a)  and  substituted  for  in  Equation  (2.2b),  thereby  diagonalizing 
the  mass  matrix.  The  ratio  Mq/M  is  denoted  and  the  equations  of  motion  are  thus 


MX  +  CX-  CaXa  +  KX-  KaXa  =  F-TX^  +  TaX^ 

MaXa  -  a^CX  +  (1  +  a'^)CaXa  -  a^KX  +  (1  +  a^)KaXa 

=  -cc^F  -  (1  +  a2)roX3  +  a^VX^ 


(2.3a) 

(2.3b) 


Next,  the  time,  T,  is  nondimensionalized  using  the  system  natural  frequency,  ui  =  ^/K/M 
and  the  differentiation  is  now  with  respect  to  nondimensional  time,  t.  Both  Equation  (2.3a) 
and  Equation  (2.3b)  are  then  divided  by  K  and  a  similarity  transformation. 
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is  performed.  The  non-dimensional  equations  of  motion  for  the  nonlinear  mechanical  vibra¬ 
tion  absorber  are 


X  -I-  2(^x  -  2a(l  -I-  J)Ca®a  +  a(l  -t-  5)'^Xa  =  /  -  +  Q:7a(l  +  S)‘^xl 

x'a  -  2ai^x  +  2(1  +  (5)(1  -I-  a^)CaiCo  -  ax  -|-  (1  -I-  a^)(l  -I-  S)^Xa 

=  —af  -  (1  -I-  a^)(l  -I-  S)^jaxl  +  a'yx^ 

where  the  following  nondimensional  parameters  are  introduced: 


(2.5a) 

(2.5b) 


C  _  O/'t  i  _ /  .2 

M  ~  M 

^  =  2Cat<^fl 


r  2 

M  =  ^  7 

M  “ 


^  =  f 

K  J 
Ma 


(2.6) 


Tfl-  —  ,..2^2.^  Kjl  _ 


Jta  ~  7a  ^ 

Analysis  for  the  case  of  small  nonlinearity  can  be  performed  using  perturbation  techniques. 
For  practical  applications,  the  mass  of  the  absorber  is  often  constrained  such  that  a  is  a 
small  parameter.  Next,  introduce  a  small  parameter  e  as  a  bookkeeping  parameter,  and 
scale  the  coefficients  as  in  Table  2.1.  Using  this  scaling  and  including  terms  to  order  e,  the 


Table  2.1:  Scaling  for  the  Mechanical  Vibration  Absorber 


Variable 

IB 

Ca 

a  i 

<7-0 

n 

7a 

B 

Scaling 

IQ 

e 

B 

€ 

e 

non-dimensional  equations  of  motion  for  the  mechanical  vibration  absorber  are  found  to  be 

X  +  2e(^x  +  X  —  eaxa  =  e/  —  ejx^  (2.7a) 

x'a  +  2eCaXa  -  eax  -f-  (1  -I-  2e<J)Xo  =  -e'jaxl  (2.7b) 

In  Equations  (2.7a)-(2.7b),  x  represents  the  displacement  of  the  mass  and  Xq  is  proportional 
to  the  relative  displacement  (or  strokelength)  of  the  absorber.  These  are  the  variables 
that  the  design  criteria  are  based  upon,  and  thus  the  formulation  is  suited  for  design. 
If  instead,  the  absolute  displacement  of  the  absorber  or  the  linear  modes  are  used  (as 
in  many  references)  the  resulting  solutions  must  be  transformed  before  application,  and 


physical  insight  is  hampered.  Additionally,  the  form  of  the  nonlinearity  is  simpler  with 
this  formulation.  In  the  next  two  sections,  the  equations  of  motion  for  the  LR-shunted 
piezoelectric  vibration  absorber  and  for  Positive  Position  feedback  control  will  be  derived 
in  an  equivalent  form. 


2.2  LR-Shunted  Piezoelectric  Vibration  Absorber 


A  modal  model  of  a  structure  containing  piezoelectric  materials  can  be  idealized  as  shown 
in  Figure  2.2.  The  base  system  consists  of  a  mass  constrained  by  a  structural  spring,  Ks, 
and  a  piezoelectric  spring,  Kp,  arranged  in  parallel.  The  displacement  of  the  structure,  X, 
is  to  be  minimized  by  a  tuned  resonant  circuit  with  charge  Q  on  the  piezoelectric  electrodes. 
The  system  thus  has  two  degrees  of  freedom  or  4  states. 

The  linear  constitutive  equations  for  piezoelectric  materials,  simplified  for  one-dimensional 
transverse  actuation,  are  [47]: 


Es 

Tn 


> 


^  —hsi 

-/l31 


Ds 

Sii 


(2.8) 


Here,  the  standard  IEEE  notation  is  used  (i.e.  E  is  electric  field,  T  is  stress,  D  is  elec¬ 
trical  displacement,  5  is  strain,  is  electrical  permittivity,  h  is  the  piezoelectric  coupling 
constant,  and  is  the  elastic  modulus).  Assuming  a  standard  patch- like  application, 
these  equations  may  be  rewritten  in  terms  of  variables  more  convenient  for  this  study.  The 
equations  for  the  piezoelectric  spring  are  thus: 


< 


V 

F 


(2.9) 


where  V  is  the  voltage  of  the  piezoelectric  electrode,  F  is  the  force  of  the  spring,  Q  is 
the  charge  flowing  into  the  patch  electrodes,  X  is  the  displacement  of  the  spring,  Cp  is 
the  capacitance  of  the  patch  under  constant  strain,  Kp  is  the  stiffness  of  the  piezoelectric 
spring  under  constant  charge,  and  H  is  the  electro-mechanical  coupling  parameter.  Note 
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Figure  2.2:  The  piezoelectric  vibration  absorber. 

that  coefficients  in  these  equations  can  be  modified  for  more  complicated  geometries,  but 
would  assume  a  similar  form. 

Placing  an  inductive-resistive  (LR)  shunt  across  the  electrodes  of  the  piezoelectric  spring, 
the  equations  of  motion  for  the  mass  in  Figure  2.2  are: 

MX  +  CX  +  K^{X,  Q)  -I-  Ks{X)  -  HQ  =  F{t)  (2.10a) 

LQ  +  RQ+-^Q-HX  =  G{X,Q)  (2.10b) 

Cp 

Here,  L  is  the  shunt  inductance;  R,  the  shunt  resistance;  M,  the  structural  mass;  and  F{t)  is 
an  external  disturbance.  The  function,  G{X,  Q),  of  the  displacement  and  charge  represents 
the  nonlinearity  in  the  shunt,  Kp{X,  Q)  represents  the  nonlinear  piezoelectric  stiffness,  and 
Ks{X)  is  the  nonlinear  structural  spring  stiffness.  The  form  of  the  nonlinearity  is  assumed 
as 

G{X,Q)=V^Q^  K^{X,Q)=K^X  Ks{X)  =  KgX +  rx^  (2.11) 

These  equations  must  be  nondimensionalized  to  facilitate  the  scaling  required  by  the  per- 


16 


turbation  analysis  of  the  next  chapter.  The  nondimensional  equations  of  motion  are: 


X  +  2(^x  +  a;  —  (1  +  5)a\q  =  f{T)  —  jx^ 
q  +  r(l  +  S)^q  +  (1  +  S)‘^q  -  (1  +  5)a\x  =  -'yaq^ 


(2.12a) 

(2.12b) 


Here,  the  nondimensional  quantities  used  in  Equations  (2. 12a, 2. 12b)  are  defined  as  follows: 


.D  -  . 


_ _ 1 


^  -  V  M  -e  -  ^ 

^2  _  f>2  _  ^31  r  _  OJe  _  1 

0^1  -  -^31  -  k+K^T^  ^ 


t  =  Tu;^ 
q  —  VLQ 


r  =  RCl>uj^ 

x  =  Vmx 


(2.13) 


Differentiation  in  Equations  (2.12a,2.12b)  is  with  respect  to  nondimensional  time,  t.  Note 
that  Equations  (2.12a,2.12b,2.13)  differ  from  those  in  the  related  literature  [20, 22-24]  since 
the  constant  charge  (shorted)  stiffness,  K^,  of  the  piezoelectric  spring  is  used  in  place  of 
the  usual  constant  voltage  stiffness,  K^,  for  the  nondimensional  equations  used  in  previous 
research  because  it  simplifies  the  analysis.  Also  note,  that  the  coupling  term  is  the  gener¬ 
alized  electro-mechanical  coupling  coefficient  which  can  be  determined  experimentally,  as  a 
modal  quantity  [22]. 

Next,  introduce  a  small  parameter  e  as  a  bookkeeping  parameter,  and  scale  the  coeffi¬ 
cients  as  in  Table  2.2.  Using  this  scaling  and  including  terms  to  order  e,  the  non-dimensional 

Table  2.2:  Scaling  for  the  LR-Shunted  Piezoelectric  Vibration  Absorber 

Variable  C  ^  cti  ^  7  7o  / 

Scaling  e  e  e  e  e  e  e 

equations  of  motion  for  the  Mechanical  Vibration  Absorber  are  found  to  be 


X  +  2e^x  +  X  —  eaiq  =  ef{t)  —  eyx^ 

q  +  2erq  -l-  (1  -t-  2eJ)g  -  eaix  =  'yaQ^ 
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(2.14a) 

(2.14b) 


These  equations  are  of  the  same  form  as  those  for  the  mechanical  vibration  absorber  pre¬ 
sented  in  Section  2.1.  The  coordinates  x  and  q  are  proportional  to  the  system  displacement 
and  shunt  charge,  respectively,  thus  maintaining  physical  significance.  In  the  next  section, 
the  equations  of  motion  for  a  single-degree-of-fi:eedom  structme  under  positive  position 
feedback  control  will  be  derived  of  this  same  form. 


2.3  Positive  Position  Feedback 


A  modal  model  of  a  structure  containing  an  actuator  can  be  idealized  as  shown  in 
Figure  2.3.  A  single  degree  of  freedom  with  mass,  M,  viscous  damping,  C,  and  stiffness,  K, 
is  driven  by  an  external  force,  F.  The  displacement,  X,  of  M  is  controlled  by  an  actuation 
force,  U.  Considering  a  nonlinear  stiffness,  the  equations  of  motion  are 

MX  +  CX  +  K{X)X  =  F  +  U  (2.15) 

The  nonlinear  stiffness  is  of  the  form 


K{X)  =  KX  +  TX^ 

Introducing  the  usual  non-dimensional  parameters 


£  =  2Ca;  f=a;2 


(2.16) 


(2.17) 


Equation  (2.15)  can  be  nondimensionalized.  For  positive  position  feedback,  U  is  defined 


U  =  Gxc  (2.18a) 

x'c  +  -t-  =  Hx  (2.18b) 

Substituting  Equation  (2.18a)  and  Equation  (2.17)  into  Equation  (2.15),  defining  t  =  wT, 
and  dividing  the  resulting  equations  by 

X  -I-  2^x  -f  a;  —  5(1  -|-  (S)^Xc  =  f{t)  —  'yx^ 
x'c  +  2(1  +  S)CcXc  +  (1  +  -  5(1  +  S)^x  =  -jcxl 
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(2.19a) 

(2.19b) 


Figure  2.3:  An  Actuated  Structure 


The  sensor  and  controller  gains,  H  and  G,  have  been  set  equal  (as  is  conventional)  and 
defined  as  (3  =  If  nonsymmetric  feedback  is  required,  a  transformation  of  the 

states,  {x,  Xc}  can  be  introduced,  bringing  the  equations  into  this  form.  Interest  in  these 
equations  is  centered  around  primary  resonance.  Scaling  of  the  equations  for  perturbation 
analysis  is  presented  in  Table  2.3  After  scaling,  the  equations  to  order  e,  where  e  is  a  small 


Table  2.3:  Scaling  for  Positive  Position  Feedback 


Variable 

c 

Cc 

9 

7c 

7 

/ 

Scaling 

e 

e 

€ 

e 

e 

e 

bookkeeping  parameter,  are: 

X  +  2e(^x  +  x  —  egxc  —  ^f{t)  —  e'yx^ 
x'c  +  2e<^Xc  +  (1  +  2e6)xc  —  egx  =  —e'ycxl 

which  axe  of  the  same  form  as  those  in  the  last  two  sections. 


(2.20a) 

(2.20b) 
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2.4  Unified  Equation 


Comparing  Equation  (2.7a)  and  Equation  (2.7b)  with  Equation  (2.12a)  and  Equa¬ 
tion  (2.12b)  or  Equation  (2.20a)  and  Equation  (2.20b)  the  nonlinear  mechanical  vibration 
absorber,  piezoelectric  vibration  absorber,  or  positive  position  feedback  control  system  can 
be  modeled  with  a  single  set  of  equations: 

x\  -I-  26^X1  +  xi  —  €01X2  =  ef{t)  —  eosxf  (2.21a) 

X2  -h  26^12  +  (1  -I-  2e^i)a;2  —  eoixi  =  eSix\  (2.21b) 

The  nondimensional  parameters  in  the  unified  equations  are  defined  in  Table  2.4.  Note  that 
in  each  case,  xi  is  proportional  to  the  displacement  of  the  system  and  X2  is  proportional  to 
the  “displacement”  of  the  absorber.  In  the  next  chapter,  the  response  of  this  equation  will 
be  determined  using  perturbation  techniques. 

Table  2.4:  Unification  of  Nonlinear  Vibration  Absorbers. 


Chapter  3 


Performance  of  Vibration 
Absorbers  for  IDOF  Systems 


Having  derived  the  common  model  for  the  various  vibration  absorbers,  their  response 
must  now  be  determined.  For  small  nonlinearities,  perturbation  techniques  are  useful  for 
determining  the  response.  Herein,  the  method  of  multiple  scales  [48]  is  used  to  transform  the 
system  of  two  second  order  equations  of  motion  into  four  first  order  modulation  equations. 
Fixed  points  of  the  modulation  equations  represent  periodic  solutions  of  the  system  equa¬ 
tions.  Once  the  modulation  equations  are  found,  a  series  of  systems  is  considered  which  are 
increasingly  more  difiicult  (i.e.  they  contain  more  nonlinear  parameters).  Since  the  modu¬ 
lation  equations  are  in  general  a  set  of  coupled  nonlinear  equations,  numerical  techniques, 
namely  pseudo-arclength  continuation  methods  are  used  to  solve  for  the  response.  Since 
no  available,  robust,  code  existed,  PSAL  ,  a  MATLAB®  toolbox  [49],  was  developed.  An 
overview,  with  references,  is  found  in  Appendix  A.  The  perturbation  solutions  are  comple¬ 
mented  by  numerical  time  simulations  to  illustrate  the  effects  of  the  various  parameters. 
Many  of  the  numerical  results  axe  also  repeated  by  analog  computer  experiments  which  are 
detailed  in  Chapter  5. 

The  organization  of  the  chapter  is  as  follows.  First,  the  modulation  equations  are 
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derived.  Next,  the  influences  of  the  various  parameters  are  examined.  First  the  linear 
parameters,  then  the  nonlinear,  are  discussed.  Finally,  conclusions  are  drawn  on  the  use  of 
these  results  for  design  of  nonlinear  vibration  absorbers. 

3.1  Derivation  of  the  Modulation  Equations 

The  general  system  equations  (2.21a)-  (2.21b)  describe  the  dynamics  of  the  various 
vibration  absorbers  presented  in  the  previous  section.  (It  can  be  noted  that  positive  position 
feedback  with  limited  authority  actuation  as  considered  in  the  previous  section  can  be 
regarded  as  an  electronic  vibration  absorber.)  Note  the  effects  of  nonlinearity,  damping, 
excitation,  and  linear  coupling  have  been  brought  together  at  0(e).  To  determine  the 
response,  the  method  of  multiple  scales  [48]  was  used.  Thus  we  express  nondimensional 
time,  t,  in  its  various  scales 

To  =  t  Ti=et  T2  =  eH  •••  (3.1) 

and  thus  by  the  chain  rule  And 

Dt  =  Dq  -t-  eD\  (3.2a) 

=Dl  +  2eDQDx  +  e^Uf  +  ■■■  (3.2b) 

where  the  operator  is  the  derivative  w.r.t.  Tj,.  Expanding  xi  and  X2  in  the  small 
parameter  e. 


=  xio  +  +  0{e‘^)  (3.3a) 

X2  =  X20  +  ex2i  +  0(e'^)  (3.3b) 

and  considering  harmonic  external  forcing  /(<)  =  /cosflTo,  the  ordinary  differential  equa¬ 
tions  (2.21a-b)  are  replaced  by  partial  (but  solvable)  differential  equations  of  various  orders 
in  e. 
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0(1)  Equations 

OqXio  +  xxQ  =0  (3.4a) 

■^0^20  +  ®20  =0  (3.4b) 

0(e)  Equations 

DqXu  +xii  =  -  2DqDiXiq  -  2(IiDqXiq  +  q;iX2o  -  a3a;io  +  /  cos  flTo  (3.4c) 

DqX2i  +  X21  =  —  2O0O1X20  —  2(^iDoX2o  +  oixio  —  (3.4d) 


Note  that  the  sign  of  the  nonlinear  terms  are  such  that  hardening  terms  are  positive.  By 
solving  the  0(1)  equations,  the  following  results: 

xio  =  Ai{Ti,  ...)e*^o  +  cc  X20  =  A2iTi,  ...)e*^'>  +  cc  (3.5) 

where  “cc”  denotes  the  complex  conjugate.  By  defining  fl  =  1  +  ecr,  the  forced  response 
near  resonance  can  be  investigated.  Substituting  these  results  into  Equations  (3.4c)-(3.4d), 
secular  terms  are  identified  leading  to  the  following  set  of  modulation  equations: 


2iA'^  =2i(^iAi  +  OC1A2  —  3Q!3.A^i4i  +  (3.6a) 

A 

2iA2  =2j^2.^2  "b  ociAi  —  'i5^A^A225\A2  (3.6b) 

The  prime  denotes  differentiation  w.r.t.  Ti.  These  equations  can  then  be  solved  by  per¬ 
forming  the  “cartesian”  substitution 

Ax  =  {px-  (3.7) 

and  separating  the  modulation  equations  into  real  and  imaginary  parts: 

ct  3 

p'l  =  -  CiPi  +  ga3(p?  +  qi)qi  (3.8a) 

q'l  =  o-pi  -b  yP2  -  Cl 91  -  +  9i)Pi  +  ^  (3.8b) 

P2  =  -  ^1)92  -  y9i  -  C2P2  +  g53(pi  +  92)92  (3.8c) 

g'  =  (cr  -  Si)p2  +  yPi  -  C292  -  g<J3(pi  +  92)P2  (3.8d) 
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The  fixed  points  of  the  resulting  four-  dimensional  set  of  nonlinear  equations  can  be  com¬ 
puted  using  continuation  methods. 


3.2  The  Linear  Response 

Before  examining  the  nonlinear  response,  the  linear  response  was  first  examined  for  two 
reasons.  First,  the  linear  solution  could  be  checked  against  the  perturbation  solution  (i.e. 
the  solution  of  the  modulation  equations)  thereby  verifying  the  perturbation  approach.  Sec¬ 
ond,  it  allows  a  complete  discussion  of  vibration  absorbers  to  be  presented  in  one  reference. 
As  noted  in  the  literature  review,  the  vibration  absorber  literature  is  spread  over  many 
sources  with  no  definitive  reference  (even  for  the  linear  case). 


3.2.1  Linear  Tuning  Parameters 


Two  design  parameters  exist  for  the  classic  vibration  absorber:  the  linear  coupling 
strength,  ai,  and  the  tuning  frequency,  Ji.  For  the  undamped  system,  the  modulation 
equations  reduce  to 


(3.9a) 

(3.9b) 


since  q\  =  q2  —  0.  Equations  (3.9a)-(3.9b)  can  be  solved  for  pi,  yielding: 

Pi 

2a{a-6,)-^ 

-aif 

P2  = - ^ 

2(t((7  -  ^l)  -  -2^ 


(3.10a) 

(3.10b) 


Checking,  if  cr  =  (5i,  pi  =  0  and  P2  =  which  agrees  with  the  linear  results  for  vibration 
absorbers  (found  for  example  in  [2])  to  first  order.  In  other  words,  if  one  uses  the  parameters 
defined  in  the  last  chapter  and  expands  the  complete  linear  solution  in  e,  the  results  are 
identical  to  0(e). 
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-  stable  --  Unstable  o  Start  +  Turn  x  Hopf  *  Br  anch 


Figure  3.1:  The  LSLA  response  as  Si  varies  from  -0.2  to  0.2  as  indicated  on  the  figure. 
System  Parameters:  =  0.001,  (2  —  0,  ai  =  0.1,  /  =  1.0,  03  =  S3  =  0. 


The  linear  system  performance  is  thus  perfect  at  a  =  Si  and  rapidly  declines  as  the 
excitation  frequency  moves  away  from  the  designed  absorber  frequency.  This  result  is 
easily  seen  in  Figure  3.1  which  was  computed  using  PSAL  for  very  light  damping  of  the 
main  system.  By  varying  ^i,  the  absorber’s  effective  range  can  be  moved  to  the  desired 
excitation  frequency.  The  response  of  the  system  can  be  eliminated  by  the  action  of  the 
absorber.  However,  since  ai  is  a  small  parameter,  large  absorber  response  is  required  to  do 
so. 

Since  ai  is  related  to  the  absorber  mass  (or  piezoelectric  coupling  or  control  gain)  a 
tradeoff  between  the  stroke  of  the  absorber  and  the  coupling,  ai,  must  be  made.  This  is 
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-Stable  -Unstable  o  Start  +Turn  x  Hopf  ’Branch 


Figure  3.2;  The  LSLA  response  for  various  values  of  ai  as  indicated  on  the  figure.  System 
Parameters:  =  0.001,  C2  =  0>  <^1  =  0>  /  =  10,  03  =  ^3  =  0. 

clearly  seen  in  Figure  3.2.  Another  benefit  of  increased  ai  is  also  apparent:  the  increase  in 
suppression  bandwidth. 

The  suppression  bandwidth  is  usually  defined  as  the  frequency  span  over  which  the 
response  of  the  system  is  less  than  the  response  of  the  system  for  a  static  load  of  equal 
magnitude.  For  linear  transfer  functions,  this  implies  the  response  lies  below  zero  dB.  As 
seen  in  Figure  3.2  and  Figure  3.1,  the  suppression  bandwidth  for  a  linear  system  increases 
with  ai ,  but  is  unaffected  by  .  Mathematically,  the  suppression  bandwidth  is  defined  by 
Pi  =  ±f.  Plugging  this  value  into  Equations  (3.9a)-(3.9b),  the  in-phase  (-f-)  and  out-of- 
phase  (-)  suppression  bandwidth  edges  can  be  found.  Substituting,  a  quadratic  equation  in 
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can  be  solved  by  expansion  in  small  parameters  to  yield: 

(3.11a) 

o? 

a  =  <5  +  -^  (3.11b) 

The  bandwidth  is  thus  af,  an  value.  The  larger  ai  is,  the  more  robust  the  system  is 
to  tuning  errors.  In  general,  the  largest  feasible  value  of  ai  will  yield  the  best  results. 
This  is  not  surprising,  as  more  control  power  is  usually  associated  with  better  performance. 
However,  in  practice,  ai  is  constrained  by  other  design  considerations  and  hence  a  method 
for  achieving  larger  bandwidths  and  lower  absorber  stroke-lengths  is  desirable.  One  method 
for  doing  so  is  the  introduction  of  damping  into  the  absorber. 

3.2.2  The  Influence  of  Damping 

The  damped  vibration  absorber  is  discussed  at  length  in  many  texts  [2, 50].  In  the  classic 
analysis,  the  system  is  treated  as  undamped  with  viscous  damping  in  the  absorber  only. 
DenHartog  [16, 50]  provided  the  first  analysis  with  damping.  By  varying  the  damping, 
the  resonant  response  can  be  reduced  at  the  expense  of  decreased  narrowband  vibration 
suppression,  as  seen  in  Figure  3.3.  Prom  this  figure,  the  effect  of  absorber  damping  is  to 
reduce  the  required  stroke  (especially  at  resonance)  with  reduced  suppression  bandwidth 
(frequency  range)  and  depth  (system  response  level) .  This  classic  tradeoff  is  well  recognized 
in  the  literature  [17].  As  the  damping  is  increased,  the  broadband  response  of  the  system 
approaches  a  minimum  while  the  narrowband  response  increases.  Thus,  it  is  possible  to 
divide  vibration  absorbers  into  two  clcisses,  the  narrowband  absorber  which  requires  small 
vibration  damping  and  the  broadband  absorber  with  larger  damping.  Given  a  fixed  linear 
coupling  a\,  narrowband  problems  are  those  in  which  the  excitation  lies  between  the  two 
resonant  peaks  defined  by  the  undamped  system: 

<5^1,2  =  2(^1  \/®i  +  <^i)  +  (3-12) 
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-Stable  -Unstable  o  Start  +Turn  x  Hopf  ‘Branch 


Figure  3.3;  The  LSLA  response  as  absorber  damping  C2  varies.  C2  =  0)  0.001,  0.01,  0.03, 
0.1,  0.3.  System  Parameters:  =  0.001,  ai  =  0.1,  Si  =  0,  f  =  1.0,  =  63  =  0. 


The  effect  of  damping  on  the  vibration  absorber  is  to  reduce  the  required  stroke-length  at 
the  cost  of  reduced  suppression  band.  If  the  damping  is  too  high,  the  suppression  band  is 
eliminated.  Such  a  device  is  still  useful  for  broadband  applications,  but,  since  the  scope  of 
this  study  is  limited  to  narrowband  disturbances,  won’t  be  discussed  further. 

A  second  somce  of  damping  is  system  damping  which  is  less  often  treated  in  the  litera¬ 
ture.  The  response  of  the  linear  absorber  as  the  damping  of  the  system  is  varied  is  plotted 
in  Figme  3.4.  Compared  with  absorber  damping,  system  damping  has  a  much  less  pro¬ 
nounced  effect  on  system  response.  Again,  the  resonant  response  is  reduced,  but  without 
the  large  penalty  in  absorber  performance.  The  depth  of  the  notch  is  reduced,  and  the 
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-Stable  -Unstable  o  Start  +Turn  x  Hopf  'Branch 


Figure  3.4:  The  LSLA  response  as  system  damping  Ci  varies.  Cl  =  Oj  0.001,  0.01,  0.03,  0.1, 
0.3.  System  Parameters:  C2  =  0.001,  ai  =  0.1,  Si  =  0,  f  =  1.0,  03  =  ^3  =  0. 


bandwidth  narrowed,  but  only  slightly.  System  damping  therefore  has  a  beneficial  effect 
of  reducing  resonant  response  with  only  slight  tradeoff  in  narrowband  suppression  perfor¬ 
mance.  Increasing  the  system  damping  can  thus  be  beneficial  without  sacrificing  absorber 
performance. 

An  absorber  design  methodology  is  still  required  which  can  expand  the  bandwidth  with¬ 
out  sacrificing  absorber  performance.  Such  a  device  was  presented  in  1952  by  Roberson: 
the  nonlinear  vibration  absorber. 
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3.3  Linear  Structure,  Nonlinear  Absorber 


The  nonlinear  absorber  is  formed  by  replacing  the  linear  spring  (or  the  electrical  shunt 
circuit)  with  a  nonlinear  element.  Cubic  nonlinearities  are  considered  herein  to  both  limit 
the  scope  and  because  they  are  most  often  discussed  in  the  literature.  Quadratic  plus  cubic 
nonlinearities  were  considered  in  the  initial  study  into  the  nonlinear  piezoelectric  vibration 
absorber  [46].  The  effect  of  the  quadratic  nonlinearity,  which  will  be  shown  in  the  next 
section,  was  to  soften  the  cubic  nonlinearity  in  a  manner  similar  to  that  reported  in  [4, 48] . 
The  effect  of  a  nonlinear  absorber  on  a  linear  system  will  be  illustrated  with  the  example 
of  [46]. 

3.3.1  Example:  Nonlinear  Piezoelectric  Vibration  Absorber 

In  this  example,  the  response  of  a  mass  supported  by  a  linear  spring  to  harmonic  input 
is  considered.  A  piezoelectric  spring  is  attached  in  parallel  so  a  piezoelectric  vibration  ab¬ 
sorber  is  formed  by  attaching  an  electrical  shunt  to  its  electrodes.  Both  quadratic  and  cubic 
nonlinearities  in  the  electrical  shunt  are  considered.  These  would  correspond  to  nonlinear 
capacitance  of  the  piezoceramic,  or  could  be  added  to  the  shunt  with  electrical  components, 
hoping  to  improve  performance.  Since  our  previous  derivation  considered  only  cubic  non- 
linearities,  the  derivation  of  the  modulation  equations  for  this  example  is  presented  first. 
Then  numerical  simulations  of  the  performance  are  presented. 

For  a  harmonic  input  of  frequency  fi,  the  equations  of  motion  are: 

q  +  S'^q  =  daix  -  r6q  +  a2q^  +  azq^  (3.13a) 

x  +  x  =  Saiq  +  f  cos{^lT)  (3.13b) 

A  “perfect”  tuning  of  the  linear  system  {5  =  1)  is  assumed  throughout.  According  to  the 
method  of  multiple  scales,  x,  q  and  T  are  expanded  in  power  series  with  e,  a  small  parameter, 
and  take  the  form: 

To  =  t  Ti=et  T2  =  eH  ■■■  (3.14a) 
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X  —  XQ  +  €Xi+e^X2  + ...  (3.14b) 

q  =  qo  +  eqi  +  e^q2  +  ...  (3.14c) 

The  coupled  system  of  ordinary  differential  equations  thus  becomes  a  system  of  partial 
differential  equations  in  the  different  time  scales  (To.Ti,...).  Absorber  performance  is 
analyzed  near  primary  resonance,  so  O  =  1  +  e^a.  At  this  point,  the  nonlinearities  must  be 
scaled.  The  piezoelectric  coupling  between  the  mechanical  deformation  and  the  electrical 
charge  is  typically  small  in  piezoelectric  vibration  absorbers,  on  the  order  of  0.1.  The 
effects  of  damping,  excitation,  nonlinear  stiffness,  and  coupling  effects  must  all  occur  at  the 
same  order.  Quadratic  nonlinearities  influence  the  equations  as  a  feedback  term  and  must 
therefore  be  of  lower  order  than  cubic  to  achieve  this  balance.  Hence  the  quadradic  term  is 
scaled  to  be  first  order  in  e  while  the  damping,  excitation,  cubic  stiffness,  and  coupling  terms 
are  second  order  in  e.  Substituting  Equation  (3.14b)  into  Equation  (3.13b)  and  including 
the  proper  scaling: 

(£>0  -b  eDi  +  ...)^(go  +  egi  +  ...)  +  (5^(go  +  +  •••)  = 

e^<5Q:i(a;o  +  +  •••)  “  €^r6{Do  +  eZ)i  +  ...){qo  +  ■■•)  "I" 

^<^2(90  +  ^qi  +  •■•)^  +  €^0:3(90  +  +  •■•)^  (3.15a) 

(Dq  +  eDi  +  ...)^(xo  +  €Xi  +  ...)  +  (xq  +  exi  +  ...)  = 

e'^Saiiqo  +  eqi  +  ...)  +  €^/cos(flT)  (3.15b) 

Expanding,  and  gathering  like  powers  of  e,  yields  the  following  set  of  equations: 

Order  0 


Dlqo  +  qo  =  0 

(3.16a) 

£>qXo  +  xo  =  0 

(3.16b) 

Order  1 

Djqi  +  qi  =  -OL2ql  -  2Do£>i9o 

(3.16c) 

DqXi  +  xi  ==  -2DqDiXq 

(3.16d) 
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Order  2 


■Do92  +  92  =  — 2£)oi)i9i  —  2Z)oD29o  ~  D\q^  —  r£)o9o  —  2a29o9i  ~  ct39o  +  Q:ia(§.16e) 

L>oa;2  +  3^2  =  -2DqDiXi  -  2D0D2X0  -  D\xo  -  D\xq  +  aiqo  +  f  cos(To  +  (rT^^.ief) 


Note  that  DiX  is  the  derivative  of  x  with  respect  to  Tj.  Solving  these  equations  yields  a 
solution  of  the  form 


9  =  9o  =  ^(Ti,T2)e'^«  a:  =  xq  =  A2(Ti,r2)e'^«  (3.17) 


where  the  coefficients  must  satisfy  the  following  differential  equation  (first  order  in  T2)  to 
eliminate  secular  terms: 


—  otiA^  — 


AjAi  +  irAi 


2iA'^  -  ai>li  - 


0 

0 


(3.18a) 

(3.18b) 


Note  i  =  \/^  and  the  over  bar  indicates  the  complex  conjugate.  Solutions  are  more  readily 
obtained  if  we  use  the  polar  form  and  thus  define 


Ai  =  aie‘^1,  71  =  crT2  -  I3i 

A2  =  a2e'^^  72  =  <^T2  -  02  (3.19) 

CXNL  =  303  -  T"2 

where  oi  and  02  are  the  amplitudes  of  the  nondimensional  shunt  charge  and  mass  vibration, 
respectively.  The  result  is  four  autonomous  differential  equations  in  four  unknowns: 


,  1  1  •  / 

Oi  =  “2^“^  2“^ 

(3.20a) 

/  1  3  1  .  N 

ai7i  =  -  gaiVLOi  -  2^2  cos(7i  -  72) 

(3.20b) 

“2  =  \f  sin(72)  -  ^02  sin(7i  -  72) 

(3.20c) 

0272  =  cra2  -  cos(72)  -  ^02  cos(7i  -  72) 

(3.20d) 

These  form  the  modulation  equations  in  polar  form.  The  cartesian  form  presented  in  the 
previous  section  is  more  suitable  for  numerical  analysis  since  it  does  not  contain  singulari¬ 
ties  if  Oj  passes  through  zero.  Steady  state  solutions  axe  obtained  by  setting  the  left  hand 
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side  of  these  equations  to  zero  and  solving  the  remaining  algebraic  equations  for  the  un¬ 
knowns.  The  amplitude  of  the  mass  vibration,  02,  may  then  be  found  in  terms  of  oi  using 
Equations  (3.20a)-(3.20b): 

o>2  =  a?  -  QijvLcrai  +  (3-21) 

The  frequency  response  equations  can  then  be  obtained  by  squaring  and  adding  Equa¬ 
tions  (3.20c)-(3.20d),  resulting  in  a  third  order  equation  in  af.  The  resulting  equation, 

-1-  craNL  (a?  -  4(T^)  Ci\  +  (4r^cr^  +  (ai  -  4o-^))  a?  -  Oi\p  =  0  (3.22) 

can  thus  be  solved  for  ai,  and  the  solution  used  in  Equation  (3.21)  to  find  02,  the  amplitude 
of  the  primary  mass  vibration. 

The  frequency  response  Equation  (3.22)  can  be  solved  using  numerical  techniques.  The 
results  duplicate  those  obtained  in  the  next  section,  so  instead  the  numerical  integration  of 
the  system  equations  is  presented.  The  system  of  equations  was  integrated  numerically  for 
a  range  of  excitation  levels  and  frequencies.  Several  sets  of  parameters  were  used.  Some  of 
these  results  are  presented  in  Figures  3.5-3.9.  The  broadening  of  the  suppression  bandwidth 
and  the  presence  of  strong  resonances  are  both  visible  in  Figures  3.5  and  3.6.  Both  the 
linear,  cubic  and  cubic-quadratic  responses  are  plotted.  The  gray  scale  indicates  the  rms 
magnitude  of  the  response.  The  rms  response  wets  calculated  using  1000  seconds  of  data 
starting  at  4000  seconds  when  the  system  had  reached  steady  state.  The  presence  of  larger 
damping  (Figure  3.6)  still  does  not  prevent  the  presence  of  strong  resonances,  even  for  small 
nonlinearities.  Only  light  damping  is  within  the  scope  of  this  work,  but  this  result  bears 
further  study  in  subsequent  investigations.  Looking  at  time  traces  in  Figures  3. 7-3. 9, the 
source  of  the  resonances  can  be  seen.  As  the  excitation  level  is  increased  from  0.1  in 
Figure  3.7  to  0.3  in  Figure  3.8,  a  bifurcation  occurs  and  the  response  becomes  quasiperiodic. 
As  the  excitation  further  increases  to  0.7  (Figure  3.9)  the  response  magnitude  appears  to 
be  chaotic  and  increases  drastically.  The  averaged  pointwise  dimension  [51]  was  found  to 
be  approximately  2.5  (Figure  3.10),  verifying  the  chaotic  nature  of  this  response.  Fifteen 


33 


Absorber  Stroke:  XI  System  Displacement:  X2 
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Figure  3.5:  Displacement  for  linear,  cubic  and  quadratic+cubic  absorbers  with  low  damping. 
These  are  plots  of  force  vs  frequency  with  the  third  parameter,  rms  displacement,  visualized 
by  shading.  Response  magnitudes  vary  from  -10  dB  (black)  to  40  dB  (white).  Responses 
were  calculated  using  ode45  in  MATLAB®  . 
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Absorber  stroke:  X1 
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System  Displacement:  X2 
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Figure  3.6:  Displacement  for  linear,  cubic  and  quadratic+cubic  absorbers  with  higher  damp¬ 
ing.  These  are  plots  of  force  vs  frequency  with  the  third  parameter, rms  displacement,  vi¬ 
sualized  by  shading.  Response  magnitudes  vary  from  -10  dB  (black)  to  40  dB  (white). 
Responses  were  calculated  using  ode45  in  MATLAB®  . 
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System  Response  to  1.025  rad/s  Forcing 
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Figure  3.7:  Numerical  simulation  of  the  piezoelectric  vibration  absorber  with  r  =  0.01, 
ai  =  0.1,  02  =  0.1,03  =  0.01,  Cl  =  1.025,  and  /  =  0.1.  Note  the  response  is  periodic  and 
the  absorber  performance  is  comparable  to  the  linear  case.  Note,  however,  the  presence  of 
higher  harmonics  in  the  frequency  response  function. 


36 


System  Response  to  1 .025  rad/s  Forcing 


Figure  3.8:  Numerical  simulation  of  the  piezoelectric  vibration  absorber  with  r  =  0.01, 
ai  =  0.1,  02  =  0.1,  as  =  0.01,f2  =  1.025,  and  /  =  0.3.  Note  the  response  is  quasiperiodic 
and  the  absorber  performance  is  maintained. 


System  Response  to  1 .025  rad/s  Forcing 


Figure  3.9:  Numerical  simulation  of  the  piezoelectric  vibration  absorber  with  r  =  0.01, 
Oil  =  O-lj  «2  =  0.1,  as  =  0.01,0  =  1.025,  and  /  =  0.7.  Note  the  response  is  chaotic  and  the 
absorber  performance  is  poor. 


thousand  seconds  of  data  were  sampled  at  thirty-two  times  the  driving  frequency.  The  first 
five  thousand  seconds  were  discarded  (leaving  55,202  samples).  One  thousand  points  were 
used  in  the  averaging  as  reference  points.  These  quasiperiodic  and  chaotic  regions  must 
therefore  he  avoided  in  designing  shunts  for  piezoelectric  vibration  absorbers  which  contain 
nonlinearities,  modeled  or  not. 

The  effect  of  small  nonlinearities  on  the  response  of  piezoelectric  vibration  absorbers 
indicates  that  while  increased  bandwidth  is  obtained  by  introducing  nonlinear  shunts,  the 
price  is  the  presence  of  resonant  response,  even  in  the  case  of  the  damped  piezoelectric  vibra¬ 
tion  absorber.  The  presence  of  periodic,  quasiperiodic,  and  chaotic  solutions  is  established. 
Further  numerical  analysis  of  the  effect  of  nonlinearity  is  required  to  better  understand 
this  behavior.  These  results  also  indicate  that  if  nonlinearity  is  present  in  a  linear  shunt 
design  and  unaccounted  for,  the  resulting  absorber  performance  could  degrade  without  the 
designer’s  knowledge. 

3.3.2  Analysis 

Motivated  by  the  bandwidth  increase  demonstrated  by  the  Nonlinear  Piezoelectric  Vi¬ 
bration  Absorber  study  presented  in  Section  3.3.1,  detailed  analysis  of  the  linear  structure 
with  a  nonlinear  absorber  was  undertaken.  Using  the  cartesian  form  of  the  modulation  equa¬ 
tions,  the  LSNA  (Linear  System,  Nonlinear  Absorber)  modulation  equations  are  obtained 
from  Equations  (3.8a)- (3. 8d): 


Pi  =  -  y  92  -  CiPi 

(3.23a) 

/  OCl  .  f 

9i  =  crpi  +  yT2  -  Ci9i  +  2 

(3.23b) 

P2  =  -  <^i)92  -  y9i  -  C2P2  +  ^^zipl  +  92)92 

(3.23c) 

q'2  =  {a-  6i)p2  -k  yPi  -  C292  -  ^^3(^2  +  92)P2 

(3.23d) 

The  performance  of  the  absorber  can  be  seen  in  Figure  3.11.  Multiple  frequency  response 
curves  are  plotted  for  varying  forcing  levels  as  indicated  in  the  caption.  When  /  =  0.01,  the 
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Figure  3.11:  The  LSNA  response  as  /  varies  for  S3  =  0.01  and  =  0.  /  =  0.01,  0.03,  0.1, 
0.3, 1.0.  System  Parameters:  =  0.001,  C2  =  0,  ai  =  0.1,  0:3  =  0 
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response  is  almost  linear  except  for  some  bend  in  the  left-hand  peak.  The  right  resonance 
bends  significantly  and  note  the  Hopf  bifurcation  which  takes  place  at  cr  =  0.085.  Within 
the  suppression  bandwidth,  the  response  of  the  system  is  acceptable. 

If  the  force  increases  to  0.03,  the  response  essentially  remains  the  same  with  larger 
bending  of  the  resonant  peaks.  However,  when  the  force  increases  to  0.1,  a  second  set  of 
Hopf  bifurcations  occur  at  <j  =  0.0019  and  a  =  0.0056  which  is  right  within  the  suppression 
region.  As  shown  in  the  last  section,  the  performance  of  the  vibration  absorber  deteriorates 
when  the  response  becomes  quasiperiodic.  Increasing  the  force  to  0.3  returns  to  a  single 
Hopf  bifurcation  (within  our  plotting  window)  at  cr  =  0.0142  and  the  suppression  region 
has  moved  with  a  zero  location  of  ct  =  0.03,  but  the  entire  suppression  region  is  quasiperi¬ 
odic.  This  agrees  with  the  solution  of  the  previous  section.  Thus  for  positive  values  of  J3, 
as  the  force  increases,  Hopf  bifurcations  appear  which  lead  to  poor  absorber  performance. 
Additionally,  the  suppression  band  moves  to  the  right  and  expands.  This  increased  sup¬ 
pression  bandwidth  is  of  little  use  especially  when  unstable.  The  coexistence  of  stable  (but 
poor  performance)  fixed  point  solutions  with  the  quasiperiodic  solutions  further  complicates 
absorber  performance. 

These  results  expand  on  those  reported  in  [34]  which  treated  the  absorber  as  a  com¬ 
bination  resonance  and  was  thus  unable  to  model  both  the  suppression  bandwidth  and 
the  resonant  behaviour  modeled  here.  Since  the  resonant  peaks  bend  into  the  suppression 
bandwidth,  the  model  must  account  for  resonant  and  suppressed  response.  This  ability  is 
a  strength  of  the  current  analysis.  The  treatment  of  the  absorber  problem  as  an  internal 
resonance  is  more  accurate  for  a  small  coupling  of  the  absorber  to  the  system.  As  seen 
in  Figure  3.11,  there  co-exists  both  suppressed  and  resonant  response  to  excitation  near 
the  absorber  tuning  frequency.  This  implies  that  the  response  of  the  system  is  dependent 
on  initial  conditions.  Moreover,  a  disturbance  to  the  system  could  cause  it  to  jump  from 
suppressed  to  resonant  behavior  implying  the  absorber  lacks  robustness  required  for  good 
design. 

Linear  tuning  of  the  nonlinear  vibration  absorber  is  examined  next.  In  Figures  3.12-3.13 
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Figure  3.12:  The  LSNA  response  as  /  varies  for  S3  =  0.01  and  (5i  =  0.05.  /  =  0.01,  0.03, 
0.1,  0.3, 1.0.  System  Parameters:  Ci  =  0.001,  C2  =  0,  qi  =  0.1,  aa  =  0 
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Figure  3.13:  The  LSNA  response  as  /  varies  for  ^3  =  0.01  and  =  —0.05.  /  =  0.001,  0.01 
0.03,  0.1,  0.3,  1.0.  System  Parameters:  =  0.001,  C2  =  0,  ai  =  0.1,  0:3  =  0 


the  linear  tuning  is  varied  with  otherwise  identical  system  parameters  as  Figure  3.11.  For 
<^1  >  0,  the  response  improves.  The  suppression  region  is  shifted  upwards  in  frequency, 
eliminating  the  robustness  issue.  Furthermore,  the  stability  of  the  system  improves  as 
although  the  Hopf  bifurcation  occurs  for  /  =  0.1,  it  is  outside  the  suppression  bandwidth. 
For  /  =  0.3,  the  unstable  region  occurs  on  the  edges  of  the  suppression  bandwidth,  but 
when  the  force  is  increased  to  /  =  1.0,  performance  again  suffers.  Negative  absorber  tuning, 
(5i  <  0,  has  similar  results  for  stability,  but  the  robustness  is  still  a  concern  since  multiple 
stable  fixed  points  are  present. 

Changing  the  sign  of  the  nonlinearity  produces  a  mirrored  response  as  plotted  in  Fig¬ 
ure  3.14  for  otherwise  identical  parameters  to  Figure  3.11.  The  system  thus  produces 
mirrored,  but  identical,  results  for  {a,<5i,<S3}  and  {— cr,— (5i,— iJa}.  For  g  >  (<)0,  S3  >  (<)0 
produces  an  increased  vibration  suppression  bandwidth.  Tuning  the  linear  absorber  fre¬ 
quency,  behaves  in  a  manner  similar  to  the  linear  system  although  the  effective  absorber 
frequency  shifts  with  increased  forcing  level,  an  important  consideration  for  nonlinear  ab¬ 
sorber  design. 

For  system  parameters  identical  except  for  damping  level  to  Figure  3.11,  numerical 
integration  of  the  system  equations  was  performed.  The  damping  was  increased  slightly  to 
improve  computation  time  (lighter  damping  results  in  longer  time  to  come  to  steady  state). 
The  equations  of  motion  were  integrated  for  3500  seconds  and  the  response  is  plotted  for  the 
next  500  seconds.  The  results  are  shown  for  multiple  forcing  levels  in  Figure  3.15.  Clearly, 
the  presence  of  Hopf  bifurcations,  and  the  corresponding  quasiperiodic  response  cause  a 
deterioration  in  absorber  performance  as  the  force  was  increased  from  0.1  to  0.3.  Note 
that  the  response  of  the  system  is  less  for  the  nonlinear  absorber  for  /  =  0.1,  but  greater 
for  /  =  0.3, 1.0.  These  results  are  similar  to  those  of  the  previous  section.  To  maintain 
performance  of  the  absorber,  the  design  force  must  be  limited.  The  cost  of  force  limitations 
must  be  offset  by  the  larger  suppression  bandwidth  evident  in  Figures  3.11-3.14  for  the 
nonlinear  absorber  to  be  a  viable  design  alternative. 

Calculation  of  the  suppression  bandwidth  for  the  undamped  LSNA  system  is  performed 
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Figure  3.14:  The  LSNA  response  as  /  varies  for  <^3  =  —0.01  and  =  0.  /  =  0.001,  0.01, 
0.03,  0.1,  0.3, 1.0.  System  Parameters:  Ci  =  0.001,  C2  =  0,  ai  =  0.1,  as  =  0 
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Figure  3.15:  The  LSNA  time  response  as  /  varies  for  S3  =  0.01  and  =  0.  /  =  0.01,  0.03, 
0.1,  0.3, 1.0.  Linear  response  (LSLA)  is  shown  dotted  on  the  figure  also.  System  Parameters: 
Cl  =  0.001,  C2  =  0,  ai  =  0.1  ,Q3  =  0 
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by  setting  Pi  +  qf  =  and  solving  the  modulation  equations  for  a  as  was  done  for  the 
linear  case  in  Section  3.2.  Two  equations  in  a  result: 


af  -  (3/4)(l  +  2a)aj(o-  -  Ji)  +  2(1  +  2(7^ S3  f  =0  (3.24a) 

af  +  (3/4)(l  -  2(r)ajia  -  Ji)  -  2(1  -  2a)® =0  (3.24b) 


for  the  in-phase  and  out-of-phase  solutions,  respectively.  Using  a  perturbation  approach, 
these  polynomial  equations  can  be  solved  for  a^: 


T+  =  Si  + 


=  Si  + 


S3f  ,  al  ,  ZSiSsf  ,  ZSU 


+  -:r  + 


at 


an 


2af 

SSiSsf 


+ 


2af 

3Sif^ 


+  0(e)' 


at 


2  2q!j  2af 
Similarly,  the  location  of  zero  system  response  can  be  found  as: 

SSsf 


+  0(6)® 


(3.25a) 

(3.25b) 


(T^  —  <Sl  + 


8a^ 


(3.26) 


Thus  like  the  linear  case  the  nonlinear  absorber  bandwidth  is  centered  around  a°  with 
magnitude 

f2  qx2  f4 


2  ^  9Sif‘ 

a,  +  —^  + 


(3.27) 


‘1 

which  is  in  agreement  with  the  numerical  results  presented  in  Figures  3.11-3.14.  Several 
authors  have  contended  that  only  softening  springs  lead  to  an  increased  absorber  bandwidth. 
Prom  Equation  (3.27),  the  analysis  shows  this  to  be  the  case  for  (5i  <  0  only.  Performance  is 
symmetric  for  the  case  of  (Ji  =  0.  For  <5i  >  0,  a  hardening  absorber  has  superior  performance 
when  first  order  nonlinear  efi'ects  and  small  coupling  ratio  ai  are  assumed.  If  the  forcing 
magnitude  is  larger,  increases  in  bandwidth  can  be  seen  for  either  case.  A  plot  of  the 
bandwidth  for  various  absorber  tunings  Si  can  be  seen  in  Figure  3.16.  The  nondimensional 
nonlinearity 


h  = 


8a® 


(3.28) 
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Figure  3.16:  The  LSNA  bandwidth  vs  S3  for  Si  as  indicated  on  the  figure.  System  Param¬ 
eters:  Cl  =  0)C2  =  =  0-1,  03  —  a2  =  S2  =  0 

is  used  in  this  plot.  This  nondimensional  group  arises  often  in  the  analysis  and  for  a  single 
nonlinearity  the  forcing  magnitude  /  can  be  considered  unity  and  the  effect  of  force  on  the 
system  is  to  change  the  effective  nonlinearity.  However,  for  multiple  nonlinear  parameters 
this  cannot  be  done  since  changes  in  /  would  result  in  simultaneous  changes  in  multiple 
parameters,  possibly  obscuring  important  relations. 

The  selection  of  absorber  design  parameters  can  be  made  using  these  results.  To  design 
an  absorber  which  has  a  given  bandwidth  at  a  frequency  cto  and  force  /,  one  first  solves 
Equation  (3.26)  for  ^i, 

(3.29) 
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and  substitutes  the  result  into  Equation  (3.27)  to  yield  the  bandwidth,  Aa: 

Aa  =  aj{l  +  — y-)  (3.30) 

The  implications  of  this  result  are  many.  First,  for  (Tq  =  0,  no  improvement  in  the  absorber 
bandwidth  can  be  achieved  while  maintaining  the  (undamped)  zero  response  at  cr  =  0. 
Second,  although  the  zero  frequency  can  theoretically  be  placed  anywhere,  care  must  be 
taken  if  it  is  placed  near  the  linear  resonances  of  the  system.  In  this  case,  the  system 
could  be  in  resonance  for  low  forcing  levels  only  to  achieve  good  performance  as  the  forcing 
increased.  Such  a  design  would  not  yield  acceptable  performance  in  most  instances.  Finally, 
limitations  on  the  achievable  bandwidth  will  arise  due  to  the  presence  of  Hopf  bifurcations 
and  their  corresponding  quasiperiodic  response. 

The  presence  of  nonlinearity  thus  has  both  benefits  and  costs.  As  seen  in  Equa¬ 
tion  (3.27),  the  suppression  bandwidth  of  the  system  can  be  increased.  This  comes  at 
the  expense  of  moving  the  location  of  zero  system  response  (perfect  absorption)  as  the  forc¬ 
ing  level  increases  and  the  possible  presence  of  Hopf  bifurcations  leading  to  quasiperiodic 
behavior  and  breakdown  of  absorption. 


3.4  Nonlinear  Structure,  Linear  Absorber 


A  second  case  to  consider  is  that  of  a  linear  absorber  on  a  nonlinear  structure.  The  as¬ 
sumption  of  a  linear  absorber  is  more  likely  to  be  accurate  for  either  the  electronic  vibration 
absorber  or  the  piezoelectric  vibration  absorber  where  electronic  components  comprise  the 
absorber.  Thus  even  for  large  displacement  of  the  structure,  the  response  of  the  absorber 
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can  be  considered  linear.  The  modulation  equations  in  this  case  reduce  to 


p'l  =  -(^Qi  -  y92  -  CiPi  +  +  qi)qi 

(3.31a) 

q[  =  api  +  yp2  -  Ci9i  -  ^oizipl  +  qj)pi  +  ^ 

(3.31b) 

P2  =  Si)q2  -  y9i  -  C2P2 

(3.31c) 

92  =  -(<T  -  5i)p2  -  ypi  -  C292 

(3.31d) 

As  in  the  last  section,  these  equations  are  solved  using  continuation  methods.  In  Figure  3.17, 

the  response  is  shown  for  a  variety  of  forcing  levels.  Contrary  to  the  nonlinearity  in  the 

absorber,  nonlinear  system  stiffness  affects  only  the  resonant  behavior,  not  the  suppression 

bandwidth.  Calculating  the  system  bandwidth  eus  in  the  previous  section  yields  the  identical 

result  to  a  linear  system: 

=  <Ji  +  ^  +  0{ef 

(3.32a) 

a- =  Sr +  0{ef 

(3.32b) 

li 

0 

b 

(3.32c) 

Robustness  is  again  a  concern  since  multiple  solutions  exist  for  some  frequencies  within  the 
absorber  suppression  region.  The  performance  of  the  absorber  for  various  forcing  levels  is 
shown  in  Figure  3.18.  The  time  response  for  /  =  {0.01,0.03,0.1,0.3, 1.0}  shows  system 
nonlinearity  does  not  affect  the  performance  of  the  absorber.  Both  the  nonlinear  and  linear 
response  are  shown,  but  overlayed  so  only  one  trace  is  visible.  However,  this  may  not  be 
the  case  from  all  initial  conditions  or  if  a  disturbance  was  to  shock  the  system,  causing  it 
to  respond  with  resonant  rather  than  suppressed  behavior. 

Thus,  the  effect  of  a  nonlinear  structure  on  the  performance  of  a  linear  absorber  is 
detrimental,  but  not  as  severe  as  for  the  nonlinear  absorber.  The  bandwidth  of  the  absorber 
is  unaffected  and,  although  Hopf  bifurcations  are  present,  they  occur  in  resonant  regions 
outside  of  the  operating  range  of  the  absorber.  However,  the  influence  of  nonlinear  systems 
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Figure  3.17:  The  NSLA  response  as  /  varies  for  03  =  0.01  and  =  0.  f  =  0.01,  0.03,  0.1 
0.3,  1.0.  System  Parameters:  =  0.001,  C2  =  0,  ai  =  0.1,  (5i  =  ^3  =  0 


on  performance  is  important  if  the  absorber  is  nonlinear,  which  will  be  demonstrated  in  the 
next  section. 


3.5  Nonlinear  Structure,  Nonlinear  Absorber 

The  final  case  considered  in  this  chapter  is  a  nonlinear  system  with  a  nonlinear  absorber. 
By  controlling  the  absorber  nonlinearity,  enhancements  in  performance  will  be  shown.  The 
modulation  equations  for  this  case  are: 


Pi  =  -(^Qi  -  y 92  -  CiPi  +  +  q\)gi 

(3.33a) 

q'l  =  crpi  +  yP2  -  Ci9i  -  ^a3(pi  +  ql)pi  +  ^ 

(3.33b) 

cti  3 

P2  =  -(^  -  <51)92  -  y91  -  C2P2  +  g(53(P2  +  92)92 

(3.33c) 

92  =  (<T  -  <5i)P2  +  yPl  -  C292  -  ^^3(P2  +  92)P2 

(3.33d) 

Fixed  points  for  these  equations  were  numerically  examined  using  continuation  methods 
and  direct  integration  as  in  the  previous  cases.  Results  for  various  forcing  levels  are  plotted 
in  Figures  3.19-3.22.  Note  that  only  three  forcing  levels  were  plotted  to  allow  a 

clearer  view  of  the  response.  The  results  for  ^303  >  0  indicate  that  the  robustness  issue  is 
more  severe  than  for  the  previous  cases.  The  resonance  curves  bend  more  severely  even  for 
low  forcing  levels.  The  symmetry  in  results  for  changing  signs  is  again  evident.  However 
when  Jaas  <  0,  the  response  is  improved  since  the  resonant  curves  bend  away  from  the 
suppression  region.  The  bandwidth  of  the  suppression  region  is  identical  to  the  LSNA 
result  Equation  (3.27).  The  system  nonlinearity  affects  only  the  resonant  peaks,  not  the 
absorber  suppression  region.  Heuristically,  if  absorption  is  occurring,  the  system  response 
is  suppressed,  minimizing  the  nonlinear  system  effects,  but  the  absorber  response  is  large. 
At  resonance,  both  system  and  absorber  response  are  large.  System  nonlinearities  are 
thus  important  in  resonant  response,  but  not  in  suppressed  response  while  the  absorber 
nonlinearities  are  significant  in  both  cases. 
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Figure  3.19:  The  NSNA  response  as  /  varies  for  az  =  0.01  and  S3  =  0.01.  /  =  0.01,  0.03, 
0.1.  System  Parameters:  =  0.001,  (2  =  0,  ai  =  0.1,  (5i  =  0 
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Figure  3.20:  The  NSNA  response  as  /  varies  for  03  =  0.01  and  S3  =  —0.01.  /  =  0.01,  0.03, 
0.1.  System  Parameters:  Cl  =  0.001,  (2  —  0,ai  =  0.1,  (5i  =  0 
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Figure  3.21;  The  NSNA  response  as  /  varies  for  as  =  —0.01  and  S3  =  0.01.  /  =  0.01,  0.03, 
0.1.  System  Parameters;  =  0.001,  C2  =  0,  ai  =  0.1,  (5i  =  0 
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Figure  3.23:  The  NSNA  response  as  /  varies  and  cr  =  0.0025  for  0:3  =  0.01  and  S3  =  0.01. 
System  Parameters:  Cl  =  0.001,  (2  =  0,  ai  =  0.1,  (5i  =  0 


59 


LOO  = }  eoo 


Figure  3.25:  The  NSNA  response  as  /  varies  and  o  =  0.0025  for  as  =  —0.01  and  (Js  =  0.01. 
System  Parameters:  =  0.001,  ^2  =  0,  ai  =  0.1,  ,5i  =  0 
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Time  responses  were  again  calculated  and  showed  similar  behavior  to  the  LSNA  results. 
Note  that  in  Figure  3.26  the  nonlinear  response  for  /  =  1.0  is  not  plotted.  This  case  diverged 
at  t  =  68  sec.,  becoming  unbounded.  Realistically,  other  nonlinear  terms  would  bound  this 
behavior.  Clearly  these  instabilities  must  be  avoided  for  successful  absorber  design. 

It  is  possible  to  achieve  significant  increases  in  absorber  bandwidth  with  acceptable 
stability  and  robustness.  For  the  softening  nonlinear  system  with  as  =  —0.01,  Ci  =  0.0, 
and  linear  coupling  fixed  as  ai=0,  a  nonlinear  absorber  was  designed  using  the  results 
of  this  chapter  for  a  frequency  of  1.05  rad/sec.  For  an  undamped  linear  absorber,  the 
suppression  bandwidth  would  be  0.01  rad/sec.  A  nonlinear  absorber  was  designed  with 
=  0.025  and  ^3  =  0.017.  The  absorber  damping  was  set  as  C2  =  0.002  in  both  cases.  The 
frequency  response  plots  for  /  =  0.2  for  both  the  linear  and  nonlinear  absorber  are  shown 
in  Figure  3.27.  A  close-up  of  the  suppression  region  is  plotted  in  Figure  3.28.  From  the 
plots,  the  suppression  bandwidth  (defined  at  -14dB)  has  more  than  doubled  from  0.009  to 
0.018  around  1.05  rad/sec  by  the  introduction  of  nonlinearity  in  the  absorber.  Note  that  the 
unstable  response  of  the  nonlinear  absorber  is  not  within  the  bandwidth  and  both  systems 
are  single-valued  within  the  suppression  bandwidth. 


3.6  Summary 

Using  the  method  of  multiple  scales,  the  unified  vibration  absorber  equation  can  be 
transformed  from  the  second  order  equations  of  motion  into  the  first  order  modulation 
equations.  These  equations  were  then  used  to  perform  parameter  studies  using  pseudo- 
arclength  continuation  with  PSAL,  a  MATLAB®  toolbox.  The  effects  of  the  conventional 
nonlinear  absorber  design  parameters  on  absorber  performance  have  been  considered. 

Three  key  absorber  performance  qualities  were  identified;  suppression  bandwidth,  sta¬ 
bility  (lack  of  quasiperiodic  response)  and  robustness  (lack  of  multiple  solutions  in  the 
absorber  bandwidth).  If  stability  is  lacking  within  the  absorber  bandwidth,  the  system  will 
undergo  large  amplitude  quasiperiodic  motions.  If  multiple  solutions  are  present,  absorber 
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Excitation  Frequency  Mistuning 


Figure  3.27:  The  NSNA  response  for  ^3  =  —0.01  and  (5i  =  0.025  (thick  lines)  and  the  NSLA 
response  for  (5$  =  0  and  =  0.05  (thin  lines)  System  Parameters:  Cl  =  0.01,  C2  =  0.002, 
ai  =  0.1,  as  =  -0.01,  /  =  0.2 
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-  Stable  -  Unstable  < 


Figure  3.28:  Closeup  of  the  NSNA  response 
the  NSLA  response  for  J3  =  0  and  =  0 
C2  =  0.002,  ai  =  0.1,  as  =  -0.01,  /  =  0.2 


Start  +  Turn  x  Hopf  *  Br  anch 


Frequency  Mistuning 


!  for  63  =  —0.01  and  =  0.025  (thick  lines)  and 
.05  (thin  lines)  System  Parameters:  (i  =  0.01, 


performance  becomes  dependent  on  initial  conditions.  An  impulse  to  the  system  could 
transfer  performance  from  suppressed  to  resonant  response.  Good  absorber  performance 
depends  on  achieving  acceptable  absorber  bandwidth  with  stability  and  robustness.  The 
stability  and  robustness  issues  have  received  only  passing  discussion  in  the  literature.  The 
effect  of  absorber  design  parameters  was  considered  using  four  systems:  the  LSLA,  LSNA, 
NSLA,  and  NSNA. 

First,  the  linear  parameters  and  (2)  were  examined.  Increasing  the  linear 

coupling  Q!i  was  determined  to  control  the  suppression  bandwidth.  As  large  a  value  of  ai 
as  possible  should  be  used.  Linear  absorber  frequency  tuning  is  accomplished  via  (5i  which 
controls  the  frequency  of  the  linear  absorber  suppression  bandwidth.  Damping  of  the  system 
Cl  has  little  effect  on  the  suppression  region,  but  reduces  the  resonant  response  significantly. 
Increased  absorber  damping,  on  the  other  hand,  reduces  and  even  eliminates  the  suppression 
bandwidth  while  simultaneously  reducing  resonant  response  of  the  system  and  the  absorber 
stroke.  It  thus  can  be  used  to  trade  suppression  bandwidth  for  strokelength  and  reduced 
resonant  response. 

The  second  system,  the  LSNA,  examined  the  effect  of  the  nonlinear  absorber  stiffness  (J3. 
Shifting  of  the  suppression  frequency  was  observed  as  the  forcing  level  of  the  system  varied. 
The  suppression  bandwidth  increases  with  increased  forcing  provided  >  0.  Otherwise, 
an  initial  decrease  followed  by  an  increase  in  the  suppression  bandwidth  was  noted.  The 
increased  performance  does,  however,  come  at  a  cost:  the  existence  of  quasiperiodic  and 
chaotic  motion  within  the  absorber  suppression  region,  affecting  the  stability  of  the  design. 
Additionally,  for  JiJa  <  0  multivalued  solutions  were  observed  within  the  suppression  region. 
Thus  although  the  suppression  bandwidth  increases,  the  stability  and  robustness  of  the 
absorber  enter  as  concerns. 

The  third  system,  the  NSLA,  demonstrated  that  nonlinear  system  stiffness  was  of  con¬ 
cern  for  linear  absorber  performance.  Although  quasiperiodic  response  was  observed,  it 
was  at  resonance,  outside  of  the  suppression  bandwidth.  Robustness  was  however  an  is¬ 
sue.  Multiple  solutions  were  possible  within  the  absorption  region  even  for  very  low  forcing 
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amplitudes.  However,  as  demonstrated  on  the  final  system,  the  NSNA,  use  of  a  nonlinear 
absorber  on  a  nonlinear  system  can  address  the  robustness  issue.  For  the  NSNA  system  the 
issue  of  stability  is  similar  to  that  of  the  LSNA. 

Performance  of  the  vibration  absorber  must  account  not  only  for  suppression  bandwidth, 
but  also  the  stability  and  robustness  of  the  design.  Limited  discussion  of  these  efi'ects  was 
found  in  the  literature.  Since  the  NSLA  can  have  robustness  issues,  the  designer  must  utilize 
a  nonlinear  absorber  to  achieve  the  more  favorable  behavior  of  the  NSNA  system  as  noted 
above.  Achieving  an  acceptable  absorber  design  will  result  in  a  tradeoff  between  acceptable 
absorber  bandwidth,  tuning  frequency  and  maiximum  disturbance  levels.  The  performance 
of  the  nonlinear  vibration  absorber  is  complex.  By  systematically  investigating  the  various 
effects,  insight  into  the  design  have  been  gained.  However,  few  real-world  systems  are  single- 
degree-of-freedom  and  thus  the  performance  of  absorbers  for  N-degree-of-freedom  systems 
was  next  considered. 
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Chapter  4 


Performance  of  Vibration 
Absorbers  for  NDOF  Systems 


Having  established  performance  criteria  and  a  design  approach  for  single-degree-of- 
freedom  systems,  N-degree-of-£reedom  systems  were  next  considered.  Motivated  by  the 
work  of  Slater  [12, 14],  an  approach  using  nonlinear  normal  modes  (NNM)  was  developed 
as  follows: 

i.  Find  the  system  equations  of  motion. 

ii.  Use  Nonlinear  Normal  Modes  to  transform  the  equations  of  motion  to  modal  coordi¬ 
nates. 

iii.  Design  the  vibration  absorber  in  modal  coordinates  for  the  two-  dimensional  modal 
manifold  being  excited.  (Note  that  it  is  assumed  that  the  mode  of  interest  is  not  in 
internal  resonance.) 

iv.  Evaluate  the  response  by  transforming  the  modal  coordinates  back  to  physical  coor¬ 
dinates  using  nonlinear  normal  mode  theory. 

To  illustrate  the  design  process,  the  example  of  Shaw  and  Pierre  [10]  is  used.  The  system. 
Figure  4.1,  consists  of  two  equal  masses  connected  with  one  nonlinear  and  two  linear  springs. 
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This  example  essentially  outlines  the  procedure  for  absorber  design  for  multi-degree-of- 
freedom  discrete  (or  discretized  continuous)  systems.  The  concepts  of  nonlinear  normal 
modes  will  be  reviewed,  then  applied  to  vibration  absorber  design  and  an  example  presented. 

4.1  Nonlinear  Normal  Modes 

As  discussed  in  Chapter  1,  two  definitions  of  Nonlinear  Normal  Modes  are  present  in 
the  literature:  one  due  to  Rosenberg  [5-7]  and  another  due  to  Shaw  and  Pierre  [10, 11]. 
For  this  work  the  definition  due  to  Shaw  and  Pierre  and  later  extended  by  Slater  [12, 14]  is 
assumed.  With  this  definition,  a  nonlinear  normal  mode  is  an  invariant  manifold  tangent  to 
the  linear  modal  plane  as  the  nonlinearity  in  the  system  decreases  to  zero.  In  the  absence 
of  internal  resonance,  each  mode  of  the  n-degree-of-freedom  system 

MijX'j  =  fi{xj,xj)  i,j  =  1,  ...,n  (4.1) 

is  a  two-dimensional  manifold  parametrized  as: 

Xi  =  Xi{u,v)  and  =  Xi  =  yi(u,  t')  (4.2) 

The  choice  of  u  and  v  is  cirbitrary.  Shaw  and  Pierre  use  the  first  physical  coordinate,  {xi,yi) 
whereas  Nayfeh  [8]  first  transforms  the  system  using  the  linear  eigenvectors  and  uses  {qk,Pk) 
where  pk  =  qk  to  find  the  nonlinear  normal  mode  tangent  to  the  linear  mode.  It  should 
be  noted  that  the  form  of  the  nonlinear  stiffness  on  the  manifold  is  dependent  on  the  choice 
of  reference  coordinates.  The  dynamics  described  remain  the  same  and  transformation 
between  definitions  is  possible.  For  the  problem  at  hand,  the  modal  coordinate  approach 
proved  simplest.  Using  standard  tensor  notation,  the  system  equations 

MjjXj  +  KijXj  -I-  €Ni{xj)  =  0  (4.3) 

undergo  a  linear  modal  transformation, 

Xi  =  ^ijqj  (4.4) 
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where  the  eigenvectors,  are  normalized  such  that 

“  ^sr  dTld  Kij^iy^gj  ~  (^*^) 

The  0  on  the  index  indicates  summation  is  not  performed  on  that  term.  The  modal 
equations  of  motion  are  then 

Qj  +  ^(j)9j  +  =  0  (4.6) 

where 


^ji^)  ~  9jrstQrQsQt 

is  assumed  for  the  form  of  the  nonlinearities  considered  for  this  study.  Using  the  method 
of  multiple  scales,  we  can  solve  for  the  nonlinear  normal  modes  as  in  [8].  Expanding  qj  as 
power  series  in  e 


Qj  —  Qjo  +  ^Qji  +  ••• 


(4.8) 


and  defining  time  as 

To  =  t,Ti  =  et, ...  (4.9) 

the  equations  may  be  ordered  in  e  to  yield 

D^gjo  +  ujfj)qjo  =  0  (4.10a) 

DoQji  +  =  -2DoDiqjo  -  gjrstQrOQsoQto  (4.10b) 

To  construct  the  nonlinear  normal  mode  which  tends  to  the  mode  as  e—  >  0,  let 

■A*  +  cc  j  =  k, 

0  j  k. 


Qjo  =  < 


(4.11) 


in  Equation  (4.10b): 


DoQji  +  ‘^fj)Qji  = 


-2DoDiqko  —  gkkkkQko  j  =  k 


(4.12) 


~gjkkkQkO 


3  ¥^k 
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Substituting  Equation  (4.11)  into  Equation  (4.12)  and  eliminating  secular  terms,  the  re¬ 
sponse  on  the  nonlinear  normal  mode  is 


Qj  =  * 

I  ^  CC  +  0(e}^  j  =  k 

_  +CC  +  0(e)2  j  ^  k 

(4.13) 

where 

3i 

D\Ak  =  - — 9kkkkA\Ak 

2u}k 

(4.14) 

can  be  solved  using  standard  perturbation  techniques.  To  place  this  result  in  a  more  common 

form,  express  (see[8]) 

qj  =  . 

\uk  j  =  k 

[rijfcgf  +  T2jkqkqk'^  j  ^k 

(4.15) 

with 

Tijk 

7tc)|  - 

(4.16) 

r2jk 

(4.17) 

This  solution  procedure  is  valid  unless  the  mode  is  in  internal  resonance.  For  the 
undamped  case  presented  here,  the  manifold  is  parametrized  by  the  single  variable 
This  will  not  be  true  for  nonconservative  systems.  However,  for  lightly  damped  and  forced 
systems,  the  response  of  the  system  to  harmonic  forcing  near  the  modal  frequency  will 
tend  to  the  undamped  nonlinear  mode  defined  above  [9, 14].  Using  this  result,  the  forced 
response  of  the  vibration  absorber  attached  to  a  forced  and  lightly  damped  structure  can 
be  analyzed. 

4.2  NDOF  absorber  design 

To  model  a  vibration  absorber  coupled  to  an  N-degree-of-ffeedom  system,  the  same 
approach  as  in  Chapter  2  can  be  followed.  The  equations  of  motion  in  physical  coordinates 
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are 

MijXj  +  KijXj  +  eNi{xj)  —  eAiXa  =  0  (4.18a) 

mx'a  +  kxa  +  eNaX^  —  eBiXi  =  0  (4.18b) 

Introducing  the  similarity  transformation  similar  to  the  previous  section, 

Xi  =  ^ijQj  (4.19a) 

Xa  =  -^Qa  (4.19b) 

y/m 

the  equations  can  be  rewritten  in  modal  coordinates: 

Qj  +  +  egjrstqrqsQt  -  ^otjqa  =  0  (4.20a) 

9a  +  +  e9a9o  “  =  0  (4.20b) 

As  in  Chapter  3  the  absorber  is  tuned  to  be  near  a  resonance  of  the  system.  Thus  Ua  =  (1  + 
Si)a)k  and  a  1:1  internal  resonance  is  formed  between  the  absorber  and  the  mode.  Note 
that  /3k  =  ctk-  Applying  the  method  of  multiple  scales,  the  equations  can  be  transformed 


into  a  series  of  partial  differential  equations.  The  0(e®)  equations  are: 

■Do9j0  +  =  0  (4.21a) 

Boqao  +  ‘^(fc)9o0  =  0  (4.21b) 

The  solution  for  the  mode  is 

* 

+  cc  j  —  k 

Qjo  =  (4.22a) 

[o  j  j^k 

q^o  =  +  cc  (4.22b) 

The  0(e)  equations  are: 

^oQji  "b  ^(j)9ji  ~  “IDfiDxqjQ  djrstQroQsoQto  "b  (^jQaO  (4.23a) 

^09ai  +  ^lk)^a0  =  -2DoDiqao  -  Qaqlo  +  “j9j0  (4.23b) 
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Substituting  Equations  (4.22a)-(4.22b)  into  Equations  (4.23a)-(4.23b), 


DoQji  +  = 


{-2DoDiqko  -  gkkkkqlo  +  OikQaO 
~9jkkkqk0  ^jQaO 


j  =  k 

3 


-^O^al  +  ~  ” 2Z)oDigaO  “  5a9a0  +  (^kQkO 


(4.24a) 

(4.24b) 


Therefore,  for  a  vibration  absorber  coupled  to  an  NDOF  system,  the  nonlinear  modal 
response  for  an  absorber  coupled  to  the  k***  linear  mode  is 


where 


Af,e^^kTo  +  ^|tyl|g3ta;fcTo  ^ 

^OjkkkAi Ak-OLj Ag  |  93jckk  a3  3iu,i.Tn  ,  „„ 


qa  =  +  cc 
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j  =  k 

3  k 


DiAk  = 


^igkkkk 

2u>k 


AlAk 


!^A 

2u>k  “ 


(4.25a) 

(4.25b) 


(4.26a) 

(4.26b) 


These  equations  are  identical  to  the  solution  for  a  single-degree-of-freedom  system  with  a 
vibration  absorber  presented  in  Chapter  3  in  the  absence  of  external  forcing  and  damping. 
The  solution  of  the  “absorbed”  mode  from  Chapter  3  can  be  used  to  find  {Ak,Aa)  which 
were  denoted  {Ai,A2)- 

Since  the  forced  and  lightly  damped  response  will  tend  towards  the  nonlinear  normal 
mode,  nonlinear  normal  mode  theory  can  be  modified  to  include  forcing  and  proportional 
damping.  First  Equations  (4.20a)-(4.20b)  axe  modified  to  include  forcing  and  damping: 


qj  +  +  u)^'^qj  +  ^gjrstqrQsQt  ^OijQa  —  IjQf 

(4.27a) 

q'a  +  2Ca9o  +  4qa  +  ^gaql  -  ^jqj  =  0 

(4.27b) 

if  +  4qf  =  0 

(4.27c) 
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Note  that  in  Equation  (4.27c)  the  forcing  has  been  written  as  a  harnaonic  variable  coupled 
to  each  mode  (linearly)  by  7j,  but  not  coupled  to  the  force.  Again  allow  uja  =  (1  + 
and  set  =  (1  +  (j)u}k  to  to  construct  the  response  of  the  mode.  Again  Pk  = 
Applying  the  method  of  multiple  scales,  the  equations  can  be  transformed  into  a  series  of 
partial  differential  equations.  The  0(e°)  equations  are: 


^oQjo  +  ^(j)QjO  =  0  (4.28a) 

Doqao  +  ^fk)QaO  =  0  (4.28b) 

^oQ/o  +  ^fk)9fo  —  0  (4.28c) 

The  solution  for  the  mode  is 


1  u  =  +  CC 

j  =  k 

(4.29a) 

9/0  =  " 

0 

j  7^  k 

9a0  -Ua  =  +  CC 

(4.29b) 

9/0  =  /o  =  +  CC 

(4.29c) 

The  0(e)  equations  are: 


Boqji  +  i^fj)qji  =  —2DQDiqjo  —  2C(j)W(j)T>o9iO  —  9jrstqroqsoqto  +  ocjqao  +  7j9/o  (4.30a) 

Dlqai  +  (^fk)^ao  =  -2DoDiqao  -  2CaW(fc)Z)o9oO  -  9aqlo  +  oijQjo  (4.30b) 

£>09/1  +  ^fk)9fi  =  -2T>o£>i9/o  (4.30c) 


Substituting  Equations  (4.29a)-(4.29b)  into  Equations  (4.30a)-(4.30c), 


£>o9/i  +  ^lj)Qji  —  < 


-2T>o£>iWo  -  ^CkDouo  —  9kkkkUo  +  akUa  +  7fc/o  j  =  k 
[  -9jkkkU^  +  ajUaO  j  9^  k 

•^09al  4"  ^(fc)9al  ~  2DQD\Ua  —  2^aDoUa  9a'^a 

£>09/1  +  i^fk)9fi  =  "•2£>o£>i/o  —  2cr/o 


(4.31a) 

(4.31b) 

(4.31c) 
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the  nonlinear  normal  mode  which  reduces  to  the  k**'  linear  mode  can  be  found.  Noting 


These  are  identical  to  the  equations  of  a  vibration  absorber  on  a  single-degree-of-freedom 
system  derived  in  Chapter  2. 

Having  solved  for  the  mode  solution,  the  l’^*’  mode  solution  is  now  derived  to  deter¬ 
mine  the  effect  of  the  vibration  absorber  on  the  other  modes  of  the  structure.  The  solution 
to  Equations  (4.27a)-(4.27c)  for  the  l’^*'  mode  is 


{ui  =  +  cc  j  =  I  ^  k 

(4.35a) 

0  j^l 


QaO  =  0 

g/0  =  0 
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(4.35b) 

(4.35c) 


The  0(e)  equations  of  motion  thus  are 


Dhji  +  = 


t 


I  —2DoDiui  —  2(iuJiDoui  -  gnnuf 


j  =  l^k 


Dlqai  +<^{kfqai  =  aiuiD^gfi  +uj^kfqfi 


=  0 


(4.36a) 

(4.36b) 


The  dynamics  of  the  nonlinear  mode  are  thus  unaffected  by  the  presence  of  the  vibration 
absorber  tuned  to  the  k*'^  mode  (in  the  absence  of  a  k-l  internal  resonance).  They  represent 
the  free  response  of  a  nonlinear  single-degree-of-freedom  system 


ill  +  2QuJiUi  +  u}fui  +  giiiiuf  =  0 


(4.37) 


For  steady  state  solutions  these  solutions  will  decay  to  zero  for  positive  damping.  The 
design  of  a  nonlinear  vibration  absorber  for  a  nonlinear  N- degree- of -freedom  system  therefore 
reduces  to  that  of  the  single-degree-of-freedom  design  presented  in  Chapter  3.  To  verify  this 
result,  a  nonlinear  two-degree-of-freedom  system  used  in  [10]  was  examined. 


4.3  Example:  2DOF  Shaw  &  Pierre  System 

The  two- mass,  three-spring  system  shown  in  Figure  4.1  was  used  in  [10]  and  later  by 
subsequent  authors  [9, 12, 14, 52]  to  illiistrate  the  concept  of  nonlinear  normal  modes.  Herein 
it  is  used  to  illustrate  the  design  procedure  developed  in  the  previous  section.  The  equations 
of  motion  of  the  2DOF  system  with  a  vibration  absorber  attached  to  mass  1  and  external 
forcing  applied  to  mass  2  are: 

MXi  +  CXi  -  CaXa  -I-  2KXx  -  KX2  +  NXf  -  KaXa  -  NaXl  =  0  (4.38a) 

Mx2  +  CX2  +  2KX2  -  KXx  =  F  (4.38b) 
Ma  (Xa  +Xx)FC^Xa  +  KaXa  +  NaX^  =  0  (4.38c) 
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Figure  4.1:  2D0F  Nonlinear  System 


Following  the  methodology  outlined  in  Chapter  2,  these  equations  can  be  nondimensional- 
ized  and  written 


x'l  +  2eca:i  +  2xi  —  X2  +  enxf  -  eaXa  —  euaxl  =  0  (4.39a) 

X2  +  2€CX2  +  2x2-  xi  =  f  (4.39b) 

x'a  +  2ecaXa  +  (1  +  2e<Ji)xo  +  €nax\  +  2eaxi  —  eaa;2  =  0  (4.39c) 


where  the  following  parameters  were  defined: 
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These  equations  are  of  the  form  assumed  in  the  previous  section.  Applying  the  modal 
transform 
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the  equations  are  transformed  into  (linear)  modal  coordinates: 


91  +  2Ci9i  +  9i  —  ^4  ^2)  =  ^ 

92  +  2-\/3^292  +  392  —  ~ 

9a  +  2Co9a  +  (1  +  2<S)9a  -  ~  ® 


(4.42a) 

(4,42b) 


(4.42c) 


Applying  the  results  of  the  previous  section,  the  nonlinear  normal  mode  of  interest  is: 


(4.43a) 

(4.43b) 

(4.43c) 

(4.44a) 

(4.44b) 


Here  u  and  Uq  represent  the  manifold  coordinates.  The  solution  to  these  equations  was 
examined  in  the  previous  chapter.  To  evaluate  the  performance  of  the  nonlinear  normal 
mode  approach,  time  integration  was  used.  Assuming  n  —  —0.04,  Ua  =  0.01,  Ci  =  0-01, 
ai  =  0.1i/(2),  Equations  (4.44a)-(4.44b)  were  integrated  numerically  using  MATLAB®  , 
and  compared  with  direct  integration  of  the  equations  of  motion  in  physical  coordinates 
and  the  0(e)  correct  modal  coordinates.  Various  harmonic  forcing  amplitudes  were  con¬ 
sidered  for  a  frequency  of  1.0  rad/s.  Results  are  plotted  in  Figures  4.2-4.13  for  both  the 
calculated  physical  coordinates,  Xi,  and  the  linear  modal  coordinates,  q,.  All  integrations 
were  performed  from  rest,  and  both  the  initial  response  and  the  steady  state  response  are 
plotted  as  indicated  on  the  figure  captions. 

Examining  Figures  4.2-4. 5,  several  observations  can  be  made.  First,  the  qualitative 
responses  of  the  figures  are  similar.  All  responses  appear  harmonic.  Second,  the  transient 
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Time  simulation  for  f=0.01  anclfreq=1  rad/s 


Figure  4.3:  Transient  response  of  the  2DOF  system  in  modal  coordinates  for  /  =  0.01  and 

(7  =  0. 
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Time  simulation  for  f=0.01  and  freq=1  rad/s 


Time  (s) 

Figure  4.4:  Steady-state  response  of  the  2DOF  system  in  physical  coordinates  for  f  =  0.01 
and  (7  =  0. 
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responses  for  all  three  cases  agree  well  for  both  qi  and  Qa,  but  the  present  nonlinear  normal 
mode  response  is  off.  This  makes  physical  sense  since  the  present  approach  considers  the 
solution  only  within  the  manifold  of  the  first  mode.  The  initial  conditions  are  not  on  that 
manifold  since  on  the  manifold  the  system  is  at  rest  only  if  /  =  0.  As  time  progresses,  the 
three  responses  deviate.  This  results  from  the  approximations  made.  Recall  that  the  slower 
time  scales  in  the  application  of  the  method  of  multiple  scales  were  not  present  in  this  order 
of  analysis.  The  transient  results  for  the  physical  coordinates  demonstrate  similar  results, 
except  the  motion  of  X2  is  in  agreement.  This  validates  the  assumption  that  the  second 
mode,  52)  does  not  participate  in  the  response  to  a  large  extent. 

Examining  the  steady  state  results  between  t=1000s  and  1200s,  the  agreement  of  the 
three  models  is  poor.  Physically,  the  reason  for  this  is  that  the  width  of  the  suppression 
bandwidth  is  an  O(e^)  quantity.  The  nonlinear  normal  mode  analysis  predicts  that  the 
suppression  bandwidth  lies  centered  near  1  rad/s.  Thus  examining  Figure  4.5,  the  NNM 
model  predicts  suppressed  response.  However  the  presence  of  the  second  mode  has  shifted 
the  suppression  bandwidth  slightly.  This  is  not  captmed  by  the  present  analysis.  The 
present  analysis  therefore  only  provides  qualitative  agreement  with  the  physical  model. 
To  achieve  quantitative  agreement  for  multi-mode  systems,  an  O(e^)  analysis  would  be 
required.  Note  that  both  the  second  mode  and  absorber  displacements  are  in  reasonable 
agreement.  This  is  because  these  responses  are  not  as  sensitive  to  the  slight  shifting  of  the 
suppression  bandwidth.  The  responses  in  physical  coordinates  shown  in  Figure  4.4  display 
similar  information. 

Increasing  the  forcing  level  to  /  =  0.1,  the  agreement  between  the  present  model  and 
the  physical  model  improves  as  shown  in  Figures  4. 6-4.9.  Again,  the  absorber  response  is  in 
excellent  agreement  while  the  system  response,  although  in  better  agreement,  shows  some 
improvement.  Recall  that  as  the  excitation  level  increases,  the  suppression  region  shifts 
higher  in  frequency  away  from  the  1  raff/s  forcing  responses  plotted.  The  sensitivity  of 
the  system  decreases  with  increased  distance  from  this  region,  so  improved  results  can  be 
expected.  Otherwise,  similar  observations  can  be  made  from  these  plots.  Note  that  the 
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Time  simulation  for  f=0.1  and  freq=1  rad/s 


Time  (s) 

Figure  4.9:  Steady-state  response  of  the  2DOF  system  in  modal  coordinates  for  /  = 
and  cr  =  0. 
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Time  simulation  for  f=0.3  and  freq=1  rad/s 


Figure  4.11:  Transient  response  of  the  2DOF  system  in  modal  coordinates  for  /  =  0.1  and 
cr  =  0. 
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Time  simulation  for  f=0.3  and  freq=1  rad/s 


Time  (s) 

Figure  4.12:  Steady-state  response  of  the  2DOF  system  in  physical  coordinates  for  /  =  0.1 
and  <7  =  0. 
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steady  state  is  periodic  in  all  three  models. 

Increasing  the  forcing  level  further  to  /  =  0.3  shows  even  better  agreement  as  plotted  in 
Figures  4.10-4.13.  The  agreement  between  the  physical  coordinates  for  both  transient  and 
modal  coordinates  is  improved.  Examining  the  steady  state  responses,  the  response  in  all 
three  cases  is  quasiperiodic.  Some  phase  difference  is  observed,  but  is  expected  due  to  the 
approximations  made  by  the  method  of  multiple  scales.  The  steady  state  magnitudes  are 
in  agreement.  Thus  qualitative  agreement  between  the  nonlinear  normal  mode  solution  and 
direct  integration  of  the  newtonian  equations  of  motion  is  achieved.  Quaiititative  agreement 
is  lacking  due  to  O(e^)  shifts  in  the  location  of  the  suppression  bandwidth  not  captured  by 
the  0(e)  approximations  used. 

4.4  Summary 

In  this  chapter,  the  equations  of  motion  for  the  N-degree-of-freedom  system  with  a 
vibration  absorber  were  shown  to  reduce  to  the  2-degree-of-freedom  absorber  model  through 
the  application  of  nonlinear  normal  mode  theory.  First  the  theory  of  nonlinear  normal  modes 
was  briefly  reviewed.  Then  the  theory  was  extended  to  allow  small  linear  coupling.  Finally, 
the  theory  of  nonlinear  normal  modes  for  internally  resonant  systems  was  extended  to  model 
the  forced  and  damped  N-degree-of-freedom  system  with  an  attached  vibration  absorber. 
This  theory  was  then  applied  to  an  example  commonly  used  in  the  nonlinear  normal  mode 
literature. 

The  example  has  verified  the  qualitative  correctness  of  the  nonlinear  mode  approxima¬ 
tion  to  the  dynamics  of  an  N-degree-of-freedom  system  with  a  vibration  absorber.  It  has 
also,  however,  highlighted  the  limitations  of  the  model  near  the  suppression  region.  These 
limitations  are  due  to  O(e^)  shifts  in  the  location  of  the  suppression  bandwidth  not  captured 
by  the  0(e)  approximations  in  the  model.  Improvement  should  be  possible  via  higher  order 
approximations . 
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Chapter  5 


Experimental  Validation 


Having  developed  the  model  in  Chapter  2  and  both  analytically  and  numerically  in¬ 
vestigated  the  performance  of  nonlinear  vibration  absorbers  in  Chapter  3,  experimental 
validation  was  undertaken.  As  noted  in  the  literature  review,  little  experimental  investi¬ 
gation  of  nonlinear  vibration  absorbers  has  been  reported.  Two  experimental  setups  were 
used.  First,  an  analog  computer  was  programmed  with  Equations  (2.21a)-(2.21b)  and  the 
effect  of  the  nonlinear  parameters  investigated.  The  analog  computer  was  then  wired  with 
the  two-degree-of-freedom  system  presented  as  an  example  in  Section  4.3.  In  the  second 
experimental  setup,  a  continuous  screw-mounted  beam  (see  Appendix  B  for  geometry)  was 
investigated  using  nonlinear  positive  position  feedback  (NPPF).  Since  only  the  first  mode 
could  be  excited  by  the  shaker  due  to  frequency  constraints,  and  no  internal  resonance 
was  identified  in  the  open-loop  structure,  the  beam  is  treated  as  a  single-degree-of-fi:eedom 
structure  and  only  the  first  mode  was  investigated.  From  the  results  of  Chapter  4,  the 
absorber  should  only  aflfect  this  mode.  Qualitative  agreement  between  the  experimental 
results  and  the  Single-  and  N-degree-of-freedom  models  is  therefore  investigated. 
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5.1  Analog  Computer  Experiments 


A  Model  7000  Compudyne  Analog  Computer  with  two  GP-10  units  was  used  to  in¬ 
vestigate  the  effect  of  the  nonlinear  parameters,  03  and  S3  on  system  response  for  both  a 
single  and  multiple(two)-degree-of-freedom  system.  Details  of  the  experimental  hardware 
used  may  be  found  in  Appendix  Chapter  B.  Two  experiments  were  conducted.  In  the 
first,  the  analog  computer  was  patched  with  Equations  (2.21a)-(2.21b)  to  investigate  the 
singe-degree-of- freedom  model.  In  the  second.  Equations  (4. 39a)- (4. 39c)  were  patched  and 
the  effects  of  the  vibration  absorber  on  both  the  first  and  second  mode  were  investigated. 

In  the  first  experiment,  the  system  was  fixed  with  a  softening  nonlinear  characteristic 
with  0:3  =  —0.01.  Four  cases  were  considered:  the  open  loop  system,  a  linear  vibration 
absorber,  a  nonlinear  softening  absorber,  and  a  nonlinear  hardening  absorber.  These  cases 
correspond  to  the  NSLA  and  NSNA  systems  modeled  in  Chapter  3.  Other  system  pa¬ 
rameters  were  fixed  as  follows:  ai  =  0.1,  =  0,  Cl  =  C2  =  0.001.  Time  scaling  was 

used,  increasing  the  characteristic  time  by  a  factor  of  100.  The  open-loop  system  natural 
frequency  is  therefore  expected  to  be  100  rad/s  or  approximately  15.9  Hz.  Two  types  of 
experiments  were  performed.  First,  swept  sine  tests  with  both  positive  and  negative  fre¬ 
quency  sweeps  were  performed  at  two  forcing  amplitudes,  /  =  0.01  and  /  =  0.1.  The  sweeps 
were  conducted  with  a  zoom  fft  bandwidth  of  0.4  Hz  and  a  step  size  of  1  Hz  from  1-10  Hz, 
0.05  Hz  from  10-20  Hz,  and  1  Hz  from  20-30  Hz.  Each  sweeps  takes  approximately  1  hour 
to  complete.  Next,  the  frequency  was  fixed  at  15.85  Hz  and  the  force  slowly  increased  to 
observe  system  response.  Time  histories  were  collected  to  evaluate  the  response.  From  the 
analytical  models,  a  Hopf  bifurcation  is  expected  with  quasiperiodic  response. 

5.1.1  Swept  Frequency  Tests 

The  swept  sine  results  are  plotted  in  Figure  5.1.  The  system  (Channel  (3))  and  absorber 
response  (Channel(5))  are  plotted  in  the  left  and  right  column,  respectively.  All  four  cases 
are  plotted.  In  the  first  row,  the  open  loop  response  is  shown  for  /  =  0.01  (the  lower 
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Frequency  (Hz.) 


Figure  5.1:  Experimental  Swept  Sine  response  of  the  single-degree-of-freedom  analog  com¬ 
puter  experiment.  The  system  (Channel  (3))  and  absorber  response  (Channel(5))  are  plot¬ 
ted  for  the  open  loop,  linear  absorber,  softening  absorber,  and  hardening  absorber  in  rows 
1-4. 
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line)  and  /  =  0.1.  Both  the  forward  and  backward  sweeps  are  plotted.  The  vertical  line 
near  12  Hz  (and  18  Hz  for  the  absorber)  results  from  the  well-known  jump  phenomenon. 
The  open-loop  absorber  consists  of  the  hardening  nonlinearity  absorber  which  is  excited  by 
the  system,  but  does  not  force  the  system  degree  of  freedom.  This  cannot  be  realized  for 
a  mechanical  or  piezoelectric  absorber,  but  is  possible  with  PPF  control.  Essentially,  the 
controller  is  excited  by  the  system,  but  the  feedback  is  turned  off.  Prom  these  plots,  the 
damping  of  the  open-loop  system  was  approximately  Ci  =  0.01.  This  increase  in  damping 
is  a  result  of  experimental  hardware  and  cannot  be  reduced.  Resonance  occurs  near  15.65 
Hz. 

In  the  second  row,  the  response  of  the  softening  system  with  linear  absorber  is  plotted. 
For  /  =  0.01,  the  response  appears  linear.  When  the  forcing  level  is  increased  to  /  =  0.1,  the 
resonant  peaks  begin  to  bend  lower  in  frequency.  This  is  expected  for  a  softening  system. 
Jump  phenomena  are  also  visible.  Resonant  response  is  begining  to  coexist  with  suppressed 
response.  Thus  the  performance  of  the  vibration  absorber  is  compromised,  as  predicted  by 
the  model  for  a  nonlinear  system  with  a  linear  absorber  (NSLA). 

Row  three  depicts  the  response  of  the  softening  system  (03  =  —0.02)  with  a  softening 
absorber  (S3  =  —0.01).  For  /  =  0.01,  the  response  is  similar  to  the  linear  system.  When 
the  sweep  is  performed  at  /  =  0.1,  the  forward  sweep  from  1  Hz  to  30  Hz  indicates  good 
performance.  The  suppression  frequency  has  shifted  negative  as  predicted  by  the  SDOF 
model  and  although  the  upper  resonance  has  shifted,  it  does  not  overlie  the  suppression 
region.  This  is  misleading  since  the  negative  sweep  from  30  Hz  to  1  Hz  indicates  that 
the  resonance  extends  from  the  upper  frequency  past  the  lower  frequency.  Although  not 
visible  on  the  plot,  it  was  observed  that  the  response  was  quasiperiodic  near  the  tip  of  the 
resonance  curve.  This  agrees  well  with  the  model  predictions. 

In  the  final  row  (row  four),  the  response  of  a  softening  system  with  a  hardening  absorber 
is  shown.  For  /  =  0.01,  the  response  again  appears  linear.  For  /  =  0.1,  the  sweeps 
indicate  acceptable  absorber  performance.  The  resonant  peaks  bend  outwards,  away  from 
the  suppression  region.  The  suppression  region  shifts  upward  in  frequency  and  slightly 
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broadens,  as  predicted  theoretically.  Jump  phenomena  and  multi-valued  solutions  exist,  but 
outside  of  the  absorber  suppression  region.  As  in  the  three  previous  cases,  the  experimental 
results  qualitatively  match  those  of  the  model  developed  in  Chapter  2  and  explored  in 
Chapter  3.  Quantitative  agreement  would  only  be  expected  to  0(e)  effects  and  as  discussed 
in  the  previous  chapter,  the  location  of  the  suppression  region  is  sensitive  to  O(e^)  effects. 

5.1.2  Swept  Excitation  Magnitude  Tests 

Having  qualitatively  verified  the  frequency  response  curves  of  the  model,  the  existence 
of  the  Hopf  bifurcations  was  next  investigated.  For  this,  the  frequency  was  fixed  at  15.85 
Hz  and  the  force  level  varied.  The  model  and  numerical  simulations  predict  that  as  the 
force  level  is  increased,  the  response  will  shift  from  periodic  to  quasiperiodic  and  absorber 
performance  will  deteriorate. 

Time  responses  for  /  =  0.01,  0.03,  0.1,  and  0.164  are  plotted  in  Figure  5.2.  As  the  force 
is  increased,  the  ratio  of  system  response  to  forcing  magnitude  decreases,  indicating  that  the 
absorption  frequency  is  shifting  upwards  with  increased  excitation  as  predicted.  The  DC 
offset  present  in  the  plots  is  a  result  of  the  analog  computer  and  could  not  be  eliminated.  It 
was  constant  in  the  test  for  various  forcing  levels.  When  the  forcing  increased  from  0.1  to 
0.164  volts,  a  change  in  response  occrured.  The  system  lost  its  periodic  nature  and  began 
to  oscillate  with  a  multi-frequency  response.  Examining  the  spectrum  of  the  response  in 
Figure  5.3,  this  response  is  qucisiperiodic. 

Further  increases  in  the  forcing  level  to  /  =  0.175,  0.19,  0.25,  and  0.3  volts  is  shown  in 
Figure  5.4.  Examining  the  responses  and  their  corresponding  spectrums  in  Figure  5.5,  the 
transition  of  the  system  response  from  a  quasiperiodic  to  a  chaotically  modulated  response 
are  presented  in  Figure  5.6.  As  the  forcing  level  is  increased  from  0.164  to  0.175  volts, 
the  character  of  the  response  changes  and  the  closed  periodic  curve  on  the  Poincare  section 
shown  in  Figure  5.7  begins  to  break  into  segments.  This  indicates  that  synchronization  of 
the  response  is  nearly  occmring.  Increasing  the  forcing  level  further,  the  response  appears 
to  become  chaotically  modulated.  The  spectrum  of  the  response  seen  in  the  lowest  plot  in 
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Figure  5.2:  Experimental  time  responses  of  the  single-degree-of-freedom  analog  computer 
experiment  for  /  =  0.01,  0.03,  0.1,  and  0.164.  The  system  and  absorber  displacements  are 
plotted  in  the  first  and  second  columns,  respectively. 
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Figure  5.3:  Experimental  spectrums  of  the  single-degree-of-freedom  analog  computer  ex¬ 
periment  for  /  =  0.01,  0.03,  0.1,  and  0.164.  The  power  spectral  density  of  the  system 
displacement  is  plotted. 
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Figure  5.4:  Experimental  time  responses  of  the  single-degree-of-£reedom  analog  computer 
experiment  for  /  =  0.175,  0.19,  0.25,  and  0.3.  The  system  and  absorber  displacements  are 
plotted  in  the  first  and  second  columns,  respectively. 
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Figure  5.5:  Experimental  spectrums  of  the  single-degree-of-freedom  analog  computer  ex¬ 
periment  for  /  =  0.175,  0.19,  0.25,  and  0.3.  The  power  spectral  density  of  the  system 
displacement  is  plotted. 
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Figure  5.7:  Experimental  Poincare  sections  of  the  single-degree-of-freedom  analog  computer 
experiment  for  all  eight  forcing  levels.  The  forcing  level  increases  down  the  first  and  then 
the  second  column,  sequentially. 
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Figure  5.5  is  broadband-its  corresponding  Poincare  section  is  space-filling.  Note  that  only 
two  of  the  four  dimensions  of  the  Poincare  section  are  plotted;  however,  given  the  nature 
of  the  time  signal,  its  spectrum,  and  its  Poincare  section,  the  response  is,  in  all  probability, 
chaotic. 

Experimental  investigation  of  a  single-degree-of-freedom  system  with  a  vibration  ab¬ 
sorber  therefore  agrees  qualitatively  with  the  model  derived  in  Chapter  2  to  0(e)  effects. 
The  quasiperiodic  response  of  the  system  and  its  transition  to  chaos  is  shown  experimentally, 
a  novel  result. 

5.1.3  Swept  Frequency  Testing  of  a  2DOF  System 

The  final  experiment  performed  on  the  analog  computer  setup  was  to  patch  the  system 
with  the  equations  of  motion  for  the  example  of  Shaw  and  Pierre  [10]  presented  in  Chapter  4. 
The  purpose  of  the  experiment  was  to  qualitatively  evaluate  the  effect  of  the  vibration 
absorber  on  the  modes  of  the  system.  The  parameters  used  were  identical  to  those  of  the 
last  chapter  with  the  exception  of  0:3  =  —0.02  instead  of -0.01.  This  corresponds  to  doubling 
the  nonlinearity,  N,  from  -0.04  to  -0.08.  This  was  done  to  facilitate  patching  three  dynamic 
systems  on  the  limited  amplifiers  of  the  analog  computer.  The  open-loop  and  closed-loop 
results  for  swept-sine  test  conducted  at  two  forcing  levels,  /  =  0.01  and  0.1,  are  shown  in 
Figme  5.8.  The  open-loop  response  is  plotted  in  row  one,  the  closed-loop  in  row  two. 
The  results  are  more  eeisily  examined  by  zooming  the  plots  near  the  two  modal  frequencies 
near  15.6  and  27.5  Hz.  Examining  the  first  mode  in  Figure  5.9,  the  effect  of  the  vibration 
absorber  is  identical  to  that  presented  in  Figure  5.1  for  the  hardening  absorber.  In  the 
second  plot.  Figure  5.10,  the  effect  of  the  vibration  absorber  tuned  to  the  first  mode  on  the 
second  mode  can  be  seen.  Perhaps  the  lack  of  effect  would  be  a  better  description.  Both 
the  open-loop  and  closed  loop  results  axe  overlaid  for  a  forcing  level  of  /  =  0.1.  These 
results  thus  qualitatively  verify  the  model  presented  in  Chapter  4.  Again  they  are  new  to 
the  literature.  Swept  forcing  tests  are  not  shown  since  the  results  duplicate  those  for  the 
single-degree-of-ffeedom  case  presented  in  the  last  section. 
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Figure  5.10:  Experimental  Swept  Sine  response  of  the  two-degree-of-freedom  analog  com¬ 
puter  experiment  near  the  second  mode.  The  system  (Channels  (3)  and  (7))  and  absorber 
(Channel(5))  responses  are  plotted  for  the  open  loop,  and  hardening  absorber  in  rows  1-2. 
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Figure  5.11:  Screw  Mounted  Beam.  See  Appendix  B  for  dimensions. 
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5.2  Screw  Mounted  Beam  Experiment 


The  second  experimental  setup  consists  of  a  beam  mounted  by  screws  on  posts  and 
actuated  by  piezoelectric  patches  surface  bonded  to  the  aluminum  beam.  Details  on  con¬ 
struction  and  hardware  used  in  the  experiment  may  be  found  in  Appendix  B.  A  diagram 
of  the  beam  is  shown  in  Figure  5.11.  Experimental  data  is  presented  first  for  the  open  loop 
(uncontrolled)  case,  then  for  two  closed  loop  cases  with  both  linear  and  nonlinear  vibration 
absorbers  implemented  as  positive  position  feedback  control  using  a  dSpace  DSP  based  real¬ 
time  controller.  Several  types  of  tests  were  run  for  a  variety  of  forcing  levels  to  investigate 
the  linear  and  nonlinear  behavior  of  the  system.  Before  PPF  control  was  implemented, 
open  loop  results  were  obtained  to  identify  the  system  under  test. 

5.2.1  Open  Loop  Results 

Three  types  of  tests  were  run  to  establish  the  baseline  model  of  the  screw-mounted  beam. 
First,  random  base  excitation  of  the  structure  was  used  to  determine  the  linear  model. 
Excitation  rms  levels  were  varied  over  several  orders  of  magnitude  to  check  for  linearity. 
The  frequency  response  transfer  functions  for  both  acceleration  and  displacement  can  be 
seen  in  Figure  5.12.  Examining  the  region  of  resonance  in  more  detail  in  Figure  5.13, 
softening  of  the  response  with  increasing  forcing  is  visible.  As  the  base  excitation  was 
increased,  resonance  tended  to  shift  lower  in  frequency. 

Since  the  random  response  showed  indications  of  nonlinear  behavior,  swept  sine  tests 
were  conducted  to  more  accurately  measure  the  system  behavior.  Frequency  steps  of  0.1 
Hz  were  taken  between  60  and  80  Hz  with  a  zoom  fit  bandwidth  of  0.4  Hz.  Four  averages  of 
the  zoomed  transfer  function  response  were  computed  at  each  point.  A  two  second  (approx. 
140  cycles)  delay  was  taken  between  frequency  steps  to  minimize  transient  effects.  Auto¬ 
ranging  was  utilized  to  ensure  accmate  data  acquisition  since  the  response  levels  varied  over 
several  orders  of  magnitude.  Each  sweep  took  approximately  1  hr.,  although  some  of  the 
nonlinear  PPF  sweeps  presented  below  took  as  long  as  4  hrs.  due  to  auto-ranging.  The 
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Sample  Rate:  512  Number  of  Averages:  32  Window:  Hanning 


Figure  5.12:  Screw  Mounted  Beam  open-loop  transfer  functions  for  various  base  acceler¬ 
ation  levels  of  0.025,  0.08,  0.25  and  0.8  grms-  Transfer  functions  are  shown  between  base 
acceleration  and  beam  acceleration  (upper  subplot)  or  beam  displacement  (lower  subplot). 


Sample  Rate:  512  Number  of  Averages:  32  Window:  Hanning 


Figure  5.13:  Screw  Mounted  Beam  open-loop  transfer  functions  for  various  base  acceleration 
levels  as  indicated  on  the  curves.  Transfer  functions  axe  shown  between  base  acceleration 
and  beam  acceleration  (upper  subplot)  or  beam  displacement  (lower  subplot). 


no 


Channels 


(2)  Beam  Accelleration 

(3)  Beam  Displacement 
(5)  PPF  Displacement 


Frequency  (Hz.)  Frequency  (Hz.)  Frequency  (Hz.) 


Figure  5.14:  Screw  Mounted  Beam  transfer  functions  for  base  acceleration  levels  of  0.01, 
0.03,  0.1  and  0.3  g’s.  Transfer  functions  between  base  acceleration  and  the  indicated  chan¬ 
nels  are  shown  for  the  open  loop,  linear  PPF  and  nonlinear  PPF  controller.  Note  that  the 
responses  are  plotted  in  volts/volt.  See  Appendix  B  for  sensor  sensitivities. 
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Figure  5.15:  Screw  Mounted  Beam  frequency  response  functions  for  various  base  accelera¬ 
tion  levels  as  indicated  on  the  cmves.  Responses  are  plotted  for  four  excitation  levels:  0.01, 
0.03,  0.1  and  0.3  g’s.  Note  that  the  responses  are  plotted  in  volts.  See  Appendix  B  for 
sensor  sensitivities. 
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open  loop  responses  are  plotted  as  transfer  functions  (Figure  5.14)  and  absolute  response 
(Figure  5.15)  in  the  first  row  of  each  figure.  Both  the  center  acceleration  and  displacement 
are  shown.  Again  similar  observations  were  made.  As  the  base  excitation  increased,  the 
resonant  response  shifts  lower  in  frequency.  This  indicates  that  a  softening  cubic  stiffness 
or  either  a  stifi'ening  or  softening  quadratic  stiffness  of  the  beam  was  present. 

Time  responses  to  a  fixed  frequency  excitation  at  various  forcing  levels  are  shown  in 
Figures  5.16-5.19.  Several  plots  are  shown  in  each  figure  for  different  forcing  levels.  First, 
portions  of  the  time  responses  on  two  channels  are  shown  in  Figure  5.16.  Both  the  system 
displacement  and  absorber  displacement  are  plotted.  For  the  open-loop  case,  the  absorber 
displacement  is  not  fed  back  to  the  piezoelectric  actuators.  Next,  the  spectrum  of  the  beam 
center  displacement  is  shown  in  Figure  5.17.  As  the  forcing  level  is  increased,  the  infiuence 
of  the  second  and  third  harmonics  varies.  The  response  is  generally  clean,  although  some 
noise  is  noted  especially  in  at  120  Hz  due  to  the  optical  displacement  sensors.  All  tests 
were  conducted  with  a  shield  to  minimize  this  noise  which  comes  from  the  fiuorescent  light 
fixtures.  X-Y  plots  for  the  time  responses  (Figure  5.16)  are  shown  in  the  first  column  of 
Figure  5.18.  In  Figure  5.19  in  the  first  column,  the  Poincare  section,  obtained  by  resampling 
the  data  at  the  forcing  frequency,  is  plotted.  These  figures  illustrate  the  noise  indicated  by 
the  spread  of  the  beam  displacement  signal  and  the  size  of  the  Poincare  section.  Instead  of 
a  single  point,  a  filled  area  is  seen,  especially  at  lower  levels.  This  is  due  to  noise  inherent 
in  experimental  work.  For  larger  forcing  levels,  the  X-Y  plot  is  an  oval,  indicating  that  the 
system  and  the  absorber  are  not  in  perfect  phase  due  to  damping  in  the  beam. 

Examination  of  the  response  spectrum  for  70.5  Hz  forcing  frequency  shows  that  both 
even  and  odd  harmonics  are  present  in  the  system  response.  However,  examination  of 
the  base  acceleration  also  showed  the  presence  of  even  and  odd  harmonics.  The  results 
were  thus  inconclusive  as  to  the  exact  nature  of  the  nonlinearity  of  the  beam  although  a 
combination  of  quadratic  and  cubic  behavior  is  suspected.  However,  as  shown  in  Chapter  3, 
the  presence  of  nonlinearity  in  the  system  only  slightly  affects  linear  absorber  performance. 
A  linear  vibration  absorber  was  next  implemented  using  positive-position-feedback  control 
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Figure  5.16:  Screw  Mounted  Beam  open-loop  time  response  in  volts  (from  top  to  bottom) 
for  x'o  =  0.005,  0.015,  0.05,  and  0.15  g.  X(l)  is  beam  displacement  and  X(2)  is  PPF 
controller  displacement.  See  Appendix  B  for  sensor  sensitivities. 


Figure  5.18:  Screw  Mounted  Beam  X-Y  time  response  in  volts  for  open-loop,  linear  and 
nonlinear  PPF  controllers,  xq  —  0.005,  0.015,  0.05,  and  0.15  g,  respectively.  X(l)  is 
beam  displacement  and  X(2)  is  PPF  controller  displacement.  See  Appendix  B  for  sensor 
sensitivities. 


to  verify  this  analytical  result. 


5.2.2  Linear  PPF 

The  PPF  control  was  implemented  using  a  dSpace  real-time  controller  as  detailed  in 
Appendix  B.  Gains  on  the  controller  were  set  using  a  trial  and  error  approach  to  maximize 
absorber  performance.  Since  the  piezoceramic  elements  on  the  beam  will  respond  over 
a  large  frequency  range,  filtering  of  the  control  signal  was  needed  to  prevent  excitation 
due  to  the  digital  signal.  The  filter  was  tuned  to  2500  Hz  which  was  half  of  the  5000  Hz 
sampling  rate.  Additionally,  an  AC  filter  was  used  on  the  beam  displacement  signal  since 
the  Philtec  optical  probes  have  a  DC  bias  in  their  output  and  played  no  role  in  the  dynamics. 
The  control  signal  filter  however  causes  a  constant  phase  shift  which  did  not  affect  the  first 
mode,  but  destabilizes  higher  modes  at  high  control  gains.  Since  the  piezoceramic  actuators 
excite  the  odd  modes,  the  third  mode  would  become  unstable  if  the  gain  is  too  high.  This 
limitation  in  gain  was  due  to  the  digital  control,  not  to  limitations  of  the  PPF  algorithm. 
To  prevent  this  problem  in  futme  work,  an  analog  PPF  filter  should  be  implemented  or  a 
piezoelectric  vibration  absorber  shunt  used  instead. 

With  the  gain  set,  linear  control  was  implemented  on  the  system.  As  with  the  open 
loop  case  three  types  of  tests  were  performed:  random,  swept  sine,  and  fixed  frequency 
excitation.  Results  from  the  random  excitation  tests  are  shown  in  Figures  5.20-5.21.  The 
response  is  slightly  different  from  the  models  of  Chapter  3,  but  generally  in  good  qualitative 
agreement.  The  asymmetry  in  the  response  is  due  to  nonlinearity  in  the  system,  agreeing 
with  the  NSLA  results  presented  in  the  Section  3.4. 

Again  steady  state  responses  to  70.5  forcing  at  0.005,  0.015,  0.05,  and  0.15  g  were 
recorded.  The  sensor  noise  is  evident  in  Figure  5.22  for  the  responses  corresponding 
to  the  lower  two  forcing  levels.  As  the  forcing  level  is  increased,  the  effect  of  noise  lessens 
due  to  increased  system  response.  The  spectrum.  Figure  5.23,  shows  some  harmonic  ef¬ 
fects,  although  less  influence  of  the  nonlinearity  is  seen  than  in  the  open  loop  response 
(Figures  5.16-5.17).  This  is  due  to  reduced  beam  displacement  as  predicted  by  the  theo- 


118 


Sample  Rate:  512  Number  of  Averages:  32  Window:  Hanning 


Frequency  (Hz.) 

Figure  5.20:  Screw  Mounted  Beam  linear  closed-loop  transfer  functions  between  base  accel¬ 
eration  and  the  indicated  channels  for  base  acceleration  levels  of  0.025,  0.08,  0.25  and  0.8 
9rms-  Channel(2)  is  beam  acceleration,  Channel  (3)  is  beam  displacement,  and  Channel  (5) 
is  PPF  filter  displacement. 
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Sample  Rate:  512  Number  of  Averages:  32  Window:  Hanning 


Frequency  (Hz.) 

Figure  5.21:  Screw  Mounted  Beam  linear  closed-loop  transfer  functions  for  base  acceleration 
levels  of  0.025,  0.08,  0.25  and  0.8  grms-  Channel(2)  is  beam  acceleration,  Channel  (3)  is 
beam  displacement,  and  Channel  (5)  is  PPF  filter  displacement. 
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3:  Screw  Mounted  Beam  with  linear  PPF  response  spectrum  for  x’o  =  0.005, 
,  and  0.15  g  (from  top  to  bottom). 


!2 


retical  results.  The  X-Y  results  and  Poincare  sections  can  be  seen  in  the  second  column 
of  Figures  5.18-5.19.  For  a  given  forcing  level,  the  system  response  is  reduced  due  to  the 
action  of  the  absorber  within  the  narrow  absorber  bandwidth.  The  linear  PPF  results  thus 
agree  with  the  single-degree-of-fireedom  model  presented  in  Chapter  3. 

5.2.3  Nonlinear  PPF 

The  experimental  effect  of  a  nonlinear  absorber  was  next  investigated.  A  nonlinear  hard¬ 
ening  term  was  added  to  the  linear  PPF  filter.  Its  level  was  set  to  maximize  the  nonlinear 
influence  without  destabilizing  higher  modes  in  the  beam.  The  reason  for  destabilization  is 
unclear  but  similar  to  increasing  the  system  gain.  The  loss  of  stability  was  for  zero  input 
force  and  it  is  hypothesized  that  the  shifting  of  phase  due  to  the  nonlinear  term  caused  the 
third  mode  to  destabilize.  Since  the  linear  gain  was  set  near  the  stability  boundary,  the 
nonlinear  phase  shift  would  cause  the  higher  mode  to  destabilize.  This  is  an  influence  of  the 
required  low-pass  filter  and  not  an  effect  of  the  PPF  vibration  absorber  per  se.  As  before, 
analog  implementation  is  recommended  to  avoid  this  problem  in  future  investigations.  If 
the  linear  gain  was  reduced,  greater  nonlinearity  was  possible.  However,  to  maintain  con¬ 
sistent  results,  this  was  not  done.  Instead,  the  nonlinear  stiffness  was  set  as  high  as  possible 
without  destabilizing  the  third  mode,  since  this  was  sufficient  to  demonstrate  the  influence 
of  nonlinearity  on  the  absorber  performance. 

Swept  sine  tests  were  conducted  and  the  results  are  shown  in  the  third  row  of  Fig¬ 
ures  5.14-5.15.  Initially,  the  nonlinear  absorber  performs  better  than  the  linear  absorber; 
system  response  is  reduced  with  less  absorber  stroke.  However,  when  the  forcing  level  is 
increased  to  0.15  g,  absorber  performance  is  degraded  and  the  system  response  is  quasiperi- 
odic.  Note  that  the  apparent  “noisy”  absorber  response  is  due  to  the  zoom  frf’s  used  in 
swept  sine  testing.  Only  the  response  around  the  driving  frequency  is  plotted,  not  the  total 
rms  of  the  system.  As  the  excitation  frequency  changes,  the  response  at  that  firequency 
varies,  giving  the  “noisy”  appearance.  Also  note  that  the  “dip”  in  system  and  absorber 
performance  on  the  0.05  g  response  curves  corresponds  to  quasiperiodic  response  at  reso- 
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nance.  Again,  the  narrowband  rms  may  decrease,  but  the  magnitude  of  system  response 
actually  increases.  The  response  simply  is  in  other  frequencies.  The  screw  beam  swept 
sine  testing  thus  agrees  with  the  analytical  models.  As  the  forcing  level  was  increased, 
quasiperiodic  system  response  occurs  due  to  a  Hopf  bifurcation  and  causes  a  breakdown  in 
absorber  performance. 

To  further  investigate  this  phenomenon,  time  responses  were  measured  at  70.5  Hz  for 
various  forcing  levels.  The  results  are  presented  in  Figure  5.24.  As  the  forcing  is  increased 
from  0.005  to  0.015  to  0.05  g’s,  absorber  performance  is  excellent.  However  when  the  force 
is  increased  to  0.15  g,  the  response  appears  quasiperiodic.  Examining  the  spectrum  in 
Figure  5.25,  the  quasiperiodic  nature  of  the  response  is  evident.  This  is  also  confirmed  by 
the  X-Y  plots  and  Poincare  sections  plotted  in  Figure  5.18  and  Figure  5.19,  respectively. 

These  results  thus  indicate  that  the  analytical  model  for  a  single-degree-of-freedom  sys¬ 
tem  with  a  vibration  absorber  is  capable  of  qualitatively  modelling  the  behavior  of  the 
screw-mounted  beam.  As  predicted,  the  nonlinearities  present  in  the  beam  did  not  signif¬ 
icantly  affect  the  performance  of  the  linear  PPF  algorithm.  However,  the  introduction  of 
absorber  nonlinearity,  while  widening  the  suppression  bandwidth,  also  introduced  Hopf  bi¬ 
furcations  resulting  in  quasiperiodic  response  and  poor  absorber  behavior.  The  nonlinearity 
in  the  absorber  is  thus  insignificant  for  low  excitation  amplitude,  beneficial  for  moderate 
excitation  amplitude,  and  detrimental  for  large  excitation  amplitude. 


5.3  Summary 

Experimental  investigation  was  conducted  on  a  single-degree-of-freedom,  a  two-degree- 
of-freedom,  and  a  continuous  system  with  vibration  absorbers  tuned  to  the  first  linear 
natural  frequency.  The  results  quantitatively  validate  the  models  presented  in  the  previ¬ 
ous  chapters.  As  the  excitation  level  increases,  experimental  verification  of  the  predicted 
quasiperiodic  and  chaotic  response  was  obtained.  Additionally,  the  need  to  tune  the  sign  of 
the  nonlinearity  of  the  absorber  opposite  to  that  of  the  system  was  observed. 


124 


Figure  5.24:  Screw  Mounted  Beam  nonlinear  PPF  time  response  in  volts  for  x'o  =  0.005 
0.015,  0.05,  and  0.15  g.  X(l)  is  beam  displacement  and  X(2)  is  PPF  controller  displacement 
See  Appendix  B  for  sensor  sensitivities. 
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Chapter  6 


Conclusions  and  Future  Work 

The  focus  of  this  work  was  to  investigate  the  performance  of  nonlinear  vibration  ab¬ 
sorbers  for  lightly  damped  multi-degree-of-fireedom  discrete  structures.  The  objectives  were 
to: 

•  Demonstrate  the  equivalence,  to  first-order  nonlinear  effects,  of  the  mechanical  vibra¬ 
tion  absorber,  the  piezoelectric  vibration  absorber,  and  positive  position  feedback. 

•  Utilize  perturbation  methods  to  find  the  general  nonlinear  response  of  the  aforemen¬ 
tioned  systems. 

•  Develop  guidelines  for  tuning  the  linear  and  nonlinear  parameters  for  vibration  ab¬ 
sorbers  for  nonlinear  single-degree-of-freedom  systems. 

•  Demonstrate  the  use  of  nonlinear  normal  modes  to  allow  modal  design  of  absorbers 
for  N-degree-of-freedom  nonlinear  structures  using  these  guidelines. 

•  Qualitatively  verify  the  models  through  experimental  investigations  for  both  single 
and  multiple  degree-of-freedom  systems. 

Each  of  these  objectives  has  been  addressed. 
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6.1  Conclusions 


In  Chapter  2,  a  general  model  for  the  mechanical  vibration  absorber,  the  Inductive- 
Resistive  Shunted  piezoelectric  vibration  absorber,  and  the  electronic  (Positive  Position 
Feedback)  vibration  absorber  was  derived.  The  model  allows  a  single  set  of  equations 
to  investigate  the  effects  of  nonlinearity  on  absorber  performance  using  a  single  dynamic 
model.  Furthermore,  this  model  was  consistent  between  systems  with  its  degrees  of  freedom 
maintaining  a  consistent  physical  meaning,  a  feature  unique  in  the  literatmre.  By  deriving 
this  model,  the  effect  of  nonlinearity  on  the  piezoelectric  and  electronic  vibration  absorber, 
an  effect  unexplored  in  the  literature  heretofore,  as  well  as  those  for  the  mechanical  vibration 
absorber  could  be  examined  using  perturbation,  numerical  and  experimental  methods  and 
the  results  applied  to  all  three  physical  implementations. 

In  Chapter  3,  the  general  model  developed  in  Chapter  2  for  the  single-degree-of-freedom 
system  with  a  vibration  absorber  were  investigated  using  both  perturbation  and  numerical 
methods.  Parameter  studies  were  performed  on  four  systems,  the  linear  system  with  a  linear 
absorber,  the  linear  system  with  a  nonlinear  absorber,  the  nonlinear  system  with  a  linear 
absorber,  and  the  nonlinear  system  with  a  nonlinear  absorber.  Three  key  performance 
qualities  were  identified;  suppression  bandwidth,  stability  and  robustness. 

The  size  of  the  suppression  bandwidth  was  found  to  depend  not  only  on  the  coupling 
parameter  (either  mass  ratio,  piezoelectric  coupling  coefficient,  or  gain)  but  also  on  the  force 
and  absorber  nonlinearity  as  a  function  of  the  parameter  group  and  the  linear  frequency 
mistuning  parameter,  <5i.  The  identification  of  these  groups  and  the  effect  of  frequency 
mistuning  is  a  novel  result  and  shows  that  for  an  absorber  tuned  above  the  open-loop 
modal  frequency,  a  hardening  absorber  is  superior  whereas  below  the  modal  frequency,  a 
softening  absorber  is  required.  The  location  in  frequency  of  the  suppression  region  was  found 
to  depend  on  the  tuning  frequency  of  the  absorber  as  well  as  the  nondimensional  nonlinear 
parameter  discussed  above.  The  combination  of  these  show  that  the  maximum  achievable 
bandwidth  increases  as  the  tuning  firequency  deviates  from  the  open-loop  modal  frequency. 
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Moreover,  it  was  determined  that  system  nonlinearity  did  not  affect  the  suppression  region 
at  the  level  of  approximation  considered  in  this  study. 

This  does  not  imply  that  system  nonlinearity  is  unimportant:  it  affects  the  second  pa¬ 
rameter,  robustness.  Robustness,  defined  in  this  study  as  the  absence  of  multiple  solutions, 
one  of  low  and  the  other  of  high  magnitude,  in  the  absorber  suppression  region.  If  the 
system  nonlinearity  is  of  the  same  type  (hardening  or  softening),  the  resonant  peaks  will 
overlie  the  suppression  region,  eliminating  system  robustness.  Additionally,  even  for  a  linear 
absorber,  robustness  is  a  concern  when  coupled  with  the  linear  mistuning.  If  the  absorber 
is  tuned  to  a  frequency  less  than  the  modal  frequency,  robustness  is  an  issue  for  hardening 
systems.  The  contrary  is  also  true.  The  end  result  is  that  for  softening  systems,  a  hardening 
nonlinear  absorber  tuned  to  a  frequency  above  the  open-loop  modal  frequency  is  desirable; 
whereas,  for  a  hardening  system,  a  softening  absorber  should  be  tuned  to  frequencies  be¬ 
low  the  modal  frequency.  This  result  clarifies  conflicting  reports  of  hardening  or  softening 
absorber  preference  in  the  literature. 

The  third  and  perhaps  most  troubling  performance  characteristic  is  stability  of  the  peri¬ 
odic  solutions  in  the  suppression  region.  As  the  excitation  level  increased,  a  Hopf  bifurcation 
was  found  to  occur,  resulting  in  quasiperiodic  response  of  the  system.  This  response  was 
determined  detrimental  to  absorber  performance  for  lightly  damped  systems.  In  the  litera¬ 
ture,  the  presence  of  higher  harmonic  behavior  has  in  some  cases  been  beneficial.  However, 
most  of  those  investigations  concerned  highly  damped  systems.  For  the  lightly  damped  sys¬ 
tems  considered  herein,  the  quasiperiodic  response  always  resulted  in  poor  performance.  For 
modal  frequencies  between  the  locations  of  linear  resonance  (i.e.  between  the  dual  natural 
frequencies  of  the  undamped  linear  vibration  absorber)  these  quasiperiodic  responses  seem 
to  limit  the  benefit  of  the  nonlinearity  to  low  forcing  levels  only.  For  higher  forcing  levels, 
increased  damping  is  required  to  suppress  this  effect,  trading  suppression  performance  for 
stability.  At  present,  no  analytic  solution  for  the  location  of  the  Hopf  bifurcations  within 
the  suppression  bandwidth  is  available  for  lightly  coupled  absorbers. 

Having  considered  the  performance  of  the  vibration  absorber  on  lightly  damped  single- 
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degree-of-freedom  systems  in  Chapter  3,  multiple-degree-of-freedom  systems  were  consid¬ 
ered  in  Chapter  4.  The  theory  of  nonlinear  normal  modes  was  extended  in  this  work  to 
include  harmonic  forcing  near  the  linear  resonance,  linear  modal  damping,  and  the  presence 
of  small  linear  coupling.  Using  these  results,  the  effect  of  vibration  absorbers  tuned  near  a 
structural  mode  was  shown  to  reduce  to  the  response  on  a  single  four  dimensional  manifold. 
The  dynamics  of  this  manifold  were  exactly  of  the  form  investigated  in  Chapters  2  and  3. 
Numerical  integration  demonstrated  some  limitations  to  the  0(e)  models  derived  in  this 
investigation.  Since  the  size  of  the  suppression  region  is  O(e^),  the  models  derived,  while 
sufficient  for  qualitative  analysis,  are  insufficient  for  quantitative  analysis.  Essentially,  since 
the  size  of  the  suppression  region  is  zero  at  the  level  of  the  analysis,  small  shifts  in  the 
location  of  this  region  can  result  from  unmodeled  effects  in  the  next  level  of  analysis.  Qual¬ 
itatively,  the  model  is  accurate  and  allows  important  insight  into  absorber  performance. 
More  importantly,  the  application  of  nonlinear  normal  mode  theory  simplifies  the  dynamics 
required  for  preliminary  analysis,  even  if  final  performance  considerations  would  require  a 
more  detailed  analysis. 

Finally,  the  accuracy  of  the  qualitative  aspects  of  the  model  were  verified  through  experi¬ 
mental  investigations  presented  in  Chapter  5.  Little  experimental  investigation  of  nonlinear 
vibration  absorbers  was  found  in  the  literature.  Through  the  analog  computer  investi¬ 
gations,  the  effects  predicted  in  Chapters  2-4  were  verified.  The  nonlinearity  shifted  the 
suppression  bandwidth  as  predicted.  The  quasiperiodic  responses  predicted  in  the  literature 
were  observed  experimentally  for  the  first  time  to  the  author’s  knowledge.  Furthermore, 
the  transition  of  the  response  to  chaotically  modulated  motion  was  observed.  The  beam 
experiment  again  verified  the  model  on  a  real-world  structure.  Quaaiperiodic  response  was 
again  observed,  although  chaotic  response  was  not,  perhaps  due  to  limitations  in  forcing 
level  and/or  feedback  gains.  Additionally,  the  effect  of  the  filtering  on  the  stability  of  higher 
modes  was  also  noted.  Finally,  the  quasiperiodic  response  was  not  a  gradual  transition,  but 
instead  an  abrupt  jump  phenomenon.  This  discovery  emphasizes  the  need  to  avoid  this 
region  since,  if  the  force  was  increased  to  the  point  where  quasiperiodic  response  occurred. 


130 


a  significant  reduction  in  the  forcing  level  was  required  to  regain  periodic  performance. 

All  of  the  goals  for  this  research  have  thus  been  accomplished.  A  consistent  model  for 
the  three  common  implementations  of  the  vibration  absorber  has  been  derived.  This  perfor¬ 
mance  of  the  nonlinear  vibration  absorber  has  been  investigated  and  qualitative  agreement 
was  found  between  the  analytical  model  and  both  numerical  and  experimental  results.  The 
performance  of  nonlinear  vibration  absorbers  for  multi-degree-of-freedom  discrete  structure 
has  been  explored  and  via  the  application  of  nonlinear  normal  modes  shown  to  reduce  to  the 
dynamics  of  the  single-degree-of-freedom  model.  Experimental  results  have  been  performed 
and  demonstrate  the  presence  of  quasiperiodic  and  chaotically-modulated  response. 


6.2  Future  Work 

This  work  would  be  incomplete  without  some  discussion  on  the  practicality  and  appli¬ 
cability  of  the  nonlinear  vibration  absorber  for  engineering  applications.  In  the  author’s 
opinion,  the  benefits  of  the  nonlinear  absorber  are  outweighed  by  the  limitations  imposed 
by  stability  concerns  for  linear  systems  at  this  time.  If  possible,  a  larger  coupling  should  be 
sought  rather  than  resorting  to  a  nonlinear  absorber.  For  a  nonlinear  structure,  the  utility 
of  the  nonlinear  absorber  may  be  to  create  a  more  robust  system  by  bending  the  resonant 
responses  away  from  the  suppression  region.  If  a  nonlinear  absorber  is  required,  the  natural 
frequency  of  the  system  should  be  shifted  so  that  the  absorber  is  operating  away  from  the 
open-loop  resonance  for  best  performance.  System  damping  should  also  be  increased  since 
the  quasiperiodic  response  can  be  delayed.  Under  these  circumstances,  it  is  possible  to 
design  a  vibration  absorber  capable  of  significantly  increasing  the  suppresion  bandwidth  for 
a  moderate  excitation  level  as  shown  in  Chapter  3.  Further  investigation  is  required  before 
nonlinear  vibration  absorbers  enter  the  vibration  engineer’s  toolbox.  Future  work  on  the 
nonlinear  vibration  absorber  can  be  divided  into  three  sections:  additional  investigations 
with  the  current  0(e)  analysis,  extension  of  this  analysis,  and  further  experimental  work. 

Using  the  current  model,  several  more  investigations  should  be  undertaken.  First,  the 
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results  of  this  study  have  been  limited  to  light  damping  since  this  case  was  least  investigated 
in  the  literature.  The  damped  nonlinear  vibration  absorber  should  be  examined.  The 
application  of  damped  vibration  absorbers  is  generally  to  broadband,  randomly  excited 
structures,  but  they  can  and  are  used  to  suppress  harmonic  excitation  as  well.  In  the 
numerical  investigation  of  the  piezoelectric  vibration  absorber  reported  in  [46],  the  presence 
of  undamped  resonant  response  in  damped  vibration  absorbers  was  seen.  This  should  be 
further  investigated.  Additionally,  more  investigation  on  the  relationship  between  damping 
and  the  minimum  excitation  level  of  the  Hopf  bifurcations  is  required  to  establish  design 
guidelines.  The  results  of  this  analysis  could  be  used  to  determine  optimal  designs  though 
the  application  of  optimization  techniques.  Such  an  analysis  should  consider  the  effects  of 
using  (static)  added  mass  or  stiffness  to  shift  system  frequencies  and  allow  larger  suppression 
bandwidths.  Any  optimization  effort  must  also  consider  the  robustness  and  stability  of  the 
final  design,  something  lacking  in  the  current  literature. 

To  allow  better  quantitative  analysis  required  to  verify  performance  and  fine-tune  de¬ 
signs,  the  0(e)  analysis  presented  herein  should  be  extended  to  the  next  level.  This  will 
allow  more  accurate  determination  of  the  suppression  region’s  location  and  width.  Addi¬ 
tionally,  the  next  level  of  analysis  would  allow  analytical  evaluation  of  preliminary  designs 
determined  using  the  current  analysis.  A  useful  result  would  be  an  analytical  determina¬ 
tion  of  the  Hopf  bifurcation  points  for  lightly  coupled  systems  such  a  more  detailed  analysis 
might  provide.  A  second  extension  of  the  analysis  would  be  to  examine  the  performance 
of  vibration  absorbers  coupled  to  a  set  of  internally  resonant  modes.  Since  periodic  and 
symmetric  structures  axe  common  in  engineering  design,  such  an  analysis  would  be  required 
to  handle  internally  resonant  modes. 

Further  experimental  investigation  of  the  application  of  nonlinear  normal  modes  to  real- 
world  structures  is  required.  The  accuracy  of  the  predicted  responses  must  be  quantitatively 
established  by  using  structures  with  analytical  models  such  as  clamped-clamped  beams.  Al¬ 
though  nonlinear,  the  screw-mounted  beam  did  not  exhibit  large  nonlinearities  at  the  forcing 
levels  obtainable  with  available  equipment.  Additionally,  although  nonlinear  Positive  Po- 
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sition  Feedback  controllers  or  inductive-resistive  piezoelectric  shunts  can  be  implemented, 
a  robust  mechanical  absorber  with  controllable  nonlinear  properties  must  be  identified  for 
those  situations  where  the  electronic  or  piezoelectric  vibration  absorber  is  impractical.  Ad¬ 
ditionally,  the  nonlinear  properties  of  piezoelectric  materials  must  to  determined  to  allow 
a  proper  nonlinear  design.  Through  the  pursuit  of  these  experimental  investigations  along 
with  the  analytical  investigations  suggested  above,  the  nonlinear  vibration  absorber  may 
be  able  to  achieve  the  level  of  practicality  required  to  join  the  engineer’s  toolchest. 
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Appendix  A 


PSAL  Overview 


PSAL  performs  pseudo-arclength  continuation  to  find  solutions  to  a  system  of  nonlinear 
equations,  {F}({x};  {A})  =  {0}.  It  is  a  robust  code  with  the  following  key  capabilities. 

•  Determination  of  fixed  points  with  or  without  time  integration  of  the  equations  {x}'  = 
{F}({x};  {A}o)  from  a  user  supplied  initial  guess  {xq}. 

•  User  function  {F}({x};{A})  is  input  along  with  its  Jacobian  in  a  single  MATLAB 
function. 

•  Continuation  is  performed  within  a  user  specified  parameter  range  with  adaptive  step- 
size  using  pseudo-arclength  continuation. 

•  Bifurcation  points  are  identified  using  efiicient  indirect  methods  and  then  solved  for 
using  direct  methods  if  desired.  Turning  Point,  Hopf,  and  Branch  points  are  sup¬ 
ported.  Branch  switching  is  not  yet  implemented. 

•  Newton-Raphson  methods  are  employed  throughout.  Second  order  differentials  are 
thus  required  for  direct  solution  of  bifurcation  points.  Numerical  computation  of  these 
derivatives  is  supported. 

This  combination  of  features  combined  with  MATLAB ’s  flexible  programing  environment 
allow  PSAL  to  serve  as  an  excellent  tool  for  parameter  study,  especially  for  small  to 
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medium  size  problems.  The  source  code  may  be  obtained  by  emailing  the  author  at 
greg . agnesOiname . com. 

The  goal  of  continuation  methods  is  to  determine  a  set  of  solutions  to  the  n  equations 
F(x;  A)  =  0  such  that  x  G  SR"  and  A  €  as  one  (or  more)  of  the  parameters,  Aj  varies  from 
an  initial  solution  {xq,  Aq}.  Some  basic  references  used  in  constructing  the  code  are  [51, 53- 
57].  The  algorithm  for  constructing  a  continuation  scheme  consists  of  three  key  portions: 
first,  a  robust  root  finder  for  systems  of  nonlinear  equations  to  find  the  initial  solution, 
{xo,  Ao};  next,  a  prediction  method  to  step  from  the  initial  solution  to  an  adjoining  solution; 
and  finally,  since  he  predicted  solution  won’t  usually  be  an  actual  point  in  the  solution  set, 
a  correction  method  to  find  the  actual  solution.  The  algorithm  continues  stepping  from 
solution  to  solution  until  a  stopping  point  is  encountered.  As  it  steps,  bifurcation  points 
must  be  identified  and  the  stability  of  the  solution  tracked.  The  algorithm  followed  must 
therefore  look  as  follows: 

i.  Find  a  starting  {k  =  0)  solution,  {xo;Ao},  to  F(x;A)  =  0. 

ii.  Evaluate  the  stability  of  this  solution. 

(a)  If  a  bifurcation  occurred  between  x^  and  Xk-i,  determine  the  type  and  use  a 
direct  method  to  find  it. 

iii.  Set  k  =  k  +  1  and  Predict  the  next  solution,  {x°;  A°} 

iv.  Correct  the  solution  with  p  iterations  and  set  Xk  =  and  A*,  =  A^. 

V.  If  {xfe,  Ajt}  is  out  of  bounds,  terminate;  otherwise  repeat  firom  step  2. 

The  theory  of  each  of  the  first  four  steps  is  briefly  explained  in  the  next  four  subsections. 
The  Fifth  step  consists  of  checking  if  A*  is  within  user  set  bounds.  For  efficiency  and  data 
structure  limitations  in  MATLAB,  the  code  actually  performs  a  user  set  block  of  iterations 
before  termination.  This  allows  a  matrix  of  fixed  size  to  be  returned  and  stored.  The  results 
are  then  stored  in  a  large  Matrix  with  an  index  identifying  the  run  (which  starting  point  it 
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began  from).  The  actual  code  does  not  terminate,  but  allows  the  user  to  begin  the  process 
again  from  a  new  starting  point.  Additionally,  plotting  commands  are  used  to  post-process 
and  plot  the  data. 

A.l  Initial  Solution 

The  initial  solution  to  the  equations  F{x;X)  =  0  is  accomplished  using  a  Newton- 
Raphson  iteration  scheme: 


II 

o 

(A.l) 

-f-  rAx^ 

(A.2) 

Fx{xi;  Xo)Ax^  =  -F{xi;Xo) 

(A.3) 

where  Fx  is  shorthand  notation  for  the  matrix  ■  The  parameter  r  is  a  relaxation 

parameter  used  to  ensure  ||jP(xfc+i, Ao)|l  <  ||F(xfc,  Ao)||.  If  not,  r  is  reduced  until  this 
relation  is  satisfied.  The  iteration  is  repeated  until 

l|F(xfe,Ao)||<£;„  (A.4) 

where  Ea  indicates  how  close  to  zero  the  norm  of  F  must  be. 

In  addition  to  the  Newton-Raphson  approach,  the  implementation  also  allows  numerical 
integration  in  time  of  the  equations  since  often  one  is  looking  for  the  steady  state  response 
to  x'  =  F{x,  A).  If  this  is  the  case,  time  integration  can  reduce  the  initial  estimate  of  the 
solution  before  Newton-Raphson  iteration  is  used,  improving  the  chance  of  convergence. 
Note  that  the  Newton-Raphson  algorithm  is  quadratically  convergent  so  it  usually  is  faster 
to  use  it  as  opposed  to  time  integration. 

A. 2  Predictor 

Once  a  solution  is  found,  a  prediction  to  the  next  fixed  point  in  {a;,A}  state-control 
space  is  needed.  It  is  tempting  to  just  increment  A  and  repeat  the  Newton-Raphson  process. 
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This  approach  will  not  yield  successful  results  near  turning  or  branch  bifurcation  points  as 
explained  in  [51,53],  for  example.  Instead  a  tangent  predictor  is  used  in  a  scheme  known 
as  arclength  continuation  [51].  Thus,  we  consider  the  arclength  s  along  the  fixed-point 
solution  curve  and  have  x  =  x{s)  and  A  =  A(s).  Differentiating  the  fixed  point  equation, 
F{x,  A)  =  0,  with  respect  to  s 

F^ix-,\)x'  +  Fxix-,X)\'  =  0  (A.5) 


a  system  of  n  equations  in  n  -I- 1  unknowns.  To  find  a  solution  we  set 

x'^x'  -k  =  As^ 


{A.6) 


which  allows  solution  for  the  predicted  step  using  an  Euler  step: 

-  ^k  +  x'As 
Afc+i  =  Xk  +  X'As 


(A.7) 

(A.8) 


This  predicted  solution  will  not  usually  be  an  exact  solution,  but  will  require  correction. 
The  size  of  the  step.  As,  is  adapted  in  the  code  using  a  method  suggested  in  [53].  The 
length  of  the  step  is  based  upon  the  number  of  corrections  needed  to  converge.  The  Euler 
step  is  slightly  modified  to 


"b  ^  '^aIIAsj; 
=  ^k  +  Ajt||Asjt 

Ask  =  Ask 

\Jdes  J 


(A.9) 

(A.IO) 

(A.ll) 


where  pk  is  the  number  of  correction  steps  needed  to  converge  to  x* ,  A* .  The  number  of 
desired  correction  steps  was  set  to  twelve,  but  can  be  adjusted  if  desired. 


A.  3  Corrector 

Once  a  prediction  A^+i}  is  determined,  it  must  be  corrected  to  satisfy  the  fixed 

point  equation.  Again  it  is  tempting  to  fix  A  and  use  the  initial  solution  iteration,  but 
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this  approach  will  face  problems  at  either  branch  or  turning  points.  To  solve  this  problem, 
pseudo-arclength  continuation  is  used.  Instead  of  seeking  a  solution  parametrized  by  A  = 
A°_(.^,  one  is  sought  which  satisfies 

g{x,  X)  =  {x-  Xkfx'k  +  (A  -  \k)X'f,  -  llxjfe,  Afe||Asfc  =  0  (A.12) 

This  is  equivalent  to  searching  for  a  solution  which  is  normal  to  the  predicted  step,  thereby 
allowing  one  to  move  around  most  turning  points.  If  small  steps  are  taken,  this  also  allows 
branch  points  to  be  included,  but  does  not  guarantee  branch  switching  will  not  occur  [53]. 
To  perform  the  correction,  an  iterative  Newton-Raphson  algorithm  is  used: 

x'k^Ax^  +  AfcAA-'  =  -g(xl^i,Xl^i) 

and 

*i+i  =  4+1  + (^15) 

4+i  =  4+i  +  ’'^v  (A.16) 

are  repeated  until  convergence.  Again,  r  is  a  relaxation  parameter.  If  the  equations  fail 
to  converge,  the  step  length.  As,  is  reduced  and  the  correction  is  repeated.  Should  this 
procedure  fail  to  converge,  a  singularity  in  the  solution  is  probable.  Prom  experience,  this 
often  happens  at  undamped  resonant  peaks.  Once  a  solution  is  found,  the  index  ,  k,  is 
updated  and  the  prediction/correction  cycle  repeated  until  termination.  Before  beginning 
the  next  prediction  cycle,  the  stability  of  the  new  fixed  point  is  investigated  to  determine 
if  bifurcations  have  occurred. 


(A.13) 

(A.14) 


A. 4  Stability  and  Bifurcation 

Two  classes  of  methods  axe  available  to  detect  and  determine  bifurcation  points:  indirect 
and  direct.  Indirect  methods  determine  the  passage  of  a  bifurcation  point  by  examining 
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data  calculated  at  the  fixed  points  determined  by  the  corrector  step.  In  general  they  are 
“cheap”  to  evaluate,  but  not  accurate  in  determining  the  location  of  the  bifurcation  point. 
They  simply  determine  that  a  bifurcation  point  was  located  between  xjt-i  and  x^-  Direct 
methods,  in  contrast,  are  more  expensive  to  evaluate  since  they  rely  on  iterative  solutions 
to  larger  systems  of  equations.  They  also  require  a  good  initial  estimate  of  the  root  location 
and  often  the  marginally  stable  eigenvector  of  the  Jacobian  (Since  bifurcation  is  marked  by 
the  passage  of  an  eigenvalue  of  the  Jacobian  into  the  right-half  plane)  for  convergence. 

PSAL  uses  a  combination  of  indirect  and  direct  methods.  Indirect  methods  are  used 
to  identify  the  presence  and  type  of  bifurcation  point  as  the  method  steps  along  the  fixed 
point  curve.  Once  a  bifurcation  point  is  identified,  its  corresponding  direct  method  is  used 
to  locate  the  point  more  exactly  if  desired.  Thus  for  quick  initial  parameter  studies  an 
approximation  can  be  used,  while  for  more  accurate  analysis,  the  direct  solution  can  be 
found.  The  drawback  to  this  approach  is  that  bifurcations  can  be  missed  if  two  bifurcation 
points  lie  between  the  steps  taken  by  the  prediction-correction  cycle.  This  possibility  can 
be  minimized  by  performing  runs  with  smaller  values  of  As,  the  prediction  step  size. 

Three  types  of  bifurcation  points  are  detected:  turning  points,  branch  points,  and  Hopf 
points.  The  first  two  are  characterized  by  having  Fx  lose  stability  through  the  origin.  Thus 
a  single  real  pole  moves  into  the  right-half  plane.  A  Hopf  bifurcation  occurs  when  two 
complex  conjugate  poles  move  into  the  right-half  plane.  Turning  points  and  branch  points 
difi'er  since  for  a  tmning  point  Fx  €  ^{Fx)  (the  null  space  of  Fx)  whereas  it  is  not  for  a 
branch  point.  By  calculating 


det  =  |Fj;|  and  DET  = 


Fx  Fx 
x'^  X' 


(A.17) 


the  presence  of  a  turning  or  branch  point  can  be  detected  since  det  changes  sign  for  turning 
points  while  DET  only  changes  at  branch  points.  Since  the  stability  of  the  fixed  point 
is  also  of  interest,  the  number  of  right-half  plane  eigenvalues  is  checked.  If  this  increases 
or  decreases  by  two,  a  Hopf  point  has  occurred.  Using  these  indirect  calculations,  the 
presence  of  bifurcation  points  is  detected.  For  quick  calculations,  the  points  are  identified 
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as  the  average  of  the  current  and  previous  fixed  points.  For  more  exact  calculations,  direct 
methods  must  be  used. 

To  directly  determine  turning  points,  an  expanded  function,  T{x,p,X)  =  0,  is  solved 
where 

>  (A.18) 

Instead  of  using  the  Newton-Raphson  method  directly,  a  faster  algorithm  from  [51]  is  used. 
The  same  algorithm  will  detect  branch  points  since  the  equations  simply  state  that  the 
Jacobian,  Fx(x,  A)  has  a  real  eigenvector  u  with  a  corresponding  eigenvalue  of  zero. 

For  Hopf  points,  H{x,p,q,\,uj)  =  0  is  solved  where 


T{x,p,X)  = 


I 


Fix-,X) 
Fx{x;  X)p 

p'^p-1 


H{x,p,q,X,oj) 


F{x-,X) 

Fx{x\X)p-uq 
<  Fx{x\X)q  +  ujp  > 
Pit  -  1 


(A.19) 


In  this  case  the  Fx  has  an  imaginary  eigenvector  p±iq  with  a  corresponding  eigenvalue  ±iuj. 
Note  the  last  two  equations  simply  normalize  the  eigenvector  such  that  the  element  is 
unity.  Again,  the  algorithm  for  solving  these  equations  is  that  in  [51]. 

Finding  the  solution  to  Equations  (A.18)-(A.19)  depends  on  an  initial  guess  for  both  the 
bifurcation  point  and  its  eigensolution.  To  determine  this,  the  eigensolution  of  the  Jacobian 
at  the  initial  guess  of  the  bifurcation  point  with  the  largest  real  portion  of  the  eigenvalue  is 
used.  Convergence  is  sometimes  still  an  issue  and  the  code  resorts  to  taking  smaller  steps 
and  using  the  indirect  scheme. 

In  summary,  the  predictor,  corrector  and  stability  algorithms  presented  form  a  robust 
MATLAB  toolbox  for  performing  a  single  parameter  study  using  pseudo-arclength  contin¬ 
uation. 
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Appendix  B 


Experimental  Setup 


This  appendix  details  the  experimental  setups  used.  Three  sections  detail  the  common 
equipment  as  well  as  the  equipment  used  for  analog  computer  simulations  and  the  screw 
mounted  beam  experiment. 

B.l  Common  Equipment 

All  experiments  were  run  with  the  general  setup  of  Figure  B.l.  The  experiment  is 
connected  to  two  I/O  systems.  One  performs  excitation  and  data  acquisition  while  the 
other  functioned  as  a  Real-time  DSP  based  controller. 

Data  acquisition  was  performed  with  two  SigLab  20-42  units  connected  to  an  Intel 
PentiumPro  200  MHz.  computer  running  Windows95.  The  SigLab  [58]  software  runs  as  a 
Matlab®  Virtual  Instruments  (Vis)  which  operate  as  an  oscilloscope,  function  generator, 
spectrum  analyzer,  network  analyzer,  or  swept  sine  analyzer.  In  addition,  long  time  record 
captiure  to  memory  or  disk  is  possible.  All  data  was  captiured  using  the  eight  differential- 
ended  inputs.  Excitation  was  via  the  four  single-ended  outputs  capable  of  producing  a 
variety  of  signals  including  random,  sine,  step,  and  arbitrary.  FFTs  are  performed  real¬ 
time  using  DSP  code  on  the  SigLab  unless  otherwise  noted.  Data  acquisition  sampling 
rate,  windowing,  and  other  relevant  parameters  are  indicated  in  the  main  text. 
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Figure  B.l;  General  Experimental  Setup. 


The  second  I/O  subsystem  is  a  DSpacea  real-time  controller.  Control  systems  are  de¬ 
signed  in  Simulink  and  downloaded  to  the  dsp  chip  for  implementation  via  the  Matlab® 
Real  Time  Workshop.  This  system  allowed  easy  implementation  of  nonlinear  real-time  con¬ 
trol,  although  its  maximum  sampling  rate  limitation  required  that  a  filter  be  used  to  smooth 
the  5  kHz  discrete  output  of  the  D/A  outputs.  Their  effect  is  discussed  in  Chapter  5. 

B.2  Analog  Computer 

For  the  analog  computer  tests  a  Compdyna  Model  7000/GP-10  system  was  used.  To 
patch  the  two-degree-of-freedom  system  with  a  vibration  absorber,  all  but  one  available 
amplifier  was  required.  Free  choice  of  the  system  parameters  was  therefore  hampered. 
Additionally  the  limitation  in  maximum  operating  speed  (time  scale)  imposes  longer  time 
requirements  on  data  acquisition.  Some  DC  offset  of  the  amplifier  channels  was  also  present, 
and  could  not  be  eliminated.  A  benefit  of  the  analog  computer  was  the  flexibility  in  adjusting 
system  parameters  in  real  time,  allowing  deeper  understanding  of  the  system  dynamics  in 
a  given  time  than  numerical  integration. 
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Figure  B.2:  Screw  Mounted  Beam.  See  text  for  dimensions. 


B.3  Screw  Mounted  Beam 

Experimental  investigations  were  conducted  on  a  screw  mounted  beam  as  seen  in  Fig¬ 
ure  B.2.  The  beam  consisted  of  an  8.5  x  1  x  0.032  inch  aluminum  substrate  with  holes 
drilled  at  1/4  in.  from  each  edge  at  all  four  corners.  Through  these  holes,  the  beam  was 
mounted  to  1/4  in.  diameter  supports  typically  used  for  spacing  electronic  circuit  boards. 
A  machine  screw  was  used  to  connect  the  beam  to  the  spacer.  The  boundary  condition 
is  not  a  classical  one,  but  of  common  engineering  practice  with  mode  shapes  similar  to  a 
clamped  boundary  condition. 

To  the  beam,  four  1.5  x  1  x  0.010  inch  nickel-electroded,  G1195  piezoceramic  sheets  were 
attached  using  MBond  2000  adhesive.  Attachment  was  performed  using  the  manufacturer’s 
instructions  for  adhering  strain  gauges  with  the  sole  modification  that  a  small  drop  of 
nickel  print  was  placed  in  the  center  of  the  sheet  to  ensure  electrical  connection  between 
the  inner  electrode  and  the  beam,  thereby  providing  a  good  ground  plane.  Polarity  of  the 
sheets  on  the  top  and  bottom  of  the  beam  were  flipped  so  a  positive  voltage  applied  to  the 
beam  resulted  in  a  bending  moment  when  the  top  sheet  expands  and  the  bottom  contracts. 
The  outer  electrodes  and  the  beam  were  then  attached  to  solder  tabs  with  fine  wire.  The 
solder  tabs  were  then  connected  to  an  amplifier  to  provide  a  control  actuator  for  the  beam. 
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Figure  B.3:  Mounting  of  the  beam. 


Maximum  amplifier  voltage  was  60  V. 

To  complete  the  beam,  a  Kistler  K-Beam  accelerometer  (Type  8302A(X)S1)  was  at¬ 
tached  at  the  center  of  the  beam  along  with  a  0.75  x  1.0  x  0.25  inch  aluminum  mass.  The 
mass  served  two  purposes.  First,  it  lowered  the  natural  frequency  of  the  system.  Second, 
the  additional  mass  causes  the  effect  of  cable  damping  from  both  the  piezoelectric  and  ac¬ 
celerometer  leads  to  be  reduced.  The  accelerometer  had  a  sensitivity  of  0.499  V/g  with  a 
maximum  of  10  g’s. 

The  beam  was  mounted  in  an  aluminum  frame  and  attached  to  an  APS  shaker  as  seen 
in  Figure  B.3.  The  aluminum  frame  had  an  additional  K-Beam  accelerometer  attached  to 
provide  a  measurement  of  the  input  motion.  The  sensitivity  of  the  accelerometer  was  499 
mv/g  with  a  maximum  of  10  g’s.  Additionally,  two  Philtec  Model  D170  displacement  probes 
were  mounted  to  provide  a  measurement  of  the  relative  displacement  of  the  beam.  One  was 
located  at  the  center  of  the  beam;  the  other  varied.  Only  the  results  from  the  center  probe 
are  reported  in  this  work.  Its  sensitivity  was  125.33  ^in/mV  with  a  maximum  of  0.1  in. 
for  linear  response.  The  center  displacement  signal  was  used  in  the  feedback  control  laws 
implemented. 
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